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1 Introduction 

in 

P\| The origins of planets, and perhaps life itself, is intrinsically linked to the 

chemistry of planet formation. In astronomy these systems are labelled as 
p ^ protoplanetary disks - disks on the incipient edge of planet formation. For our 

Sun, the rotating ball of gas and dust that collapsed to a disk has been called 
the Solar Nebula. In this chapter we will attempt to explore the chemistry of 
planet-forming disks from the perspective of knowledge gained from decades 
of solar system study, combined with our rapidly growing knowledge of extra- 
solar protoplanetary disks. This chapter is not written in the form of a review. 
+2 Rather we survey our basic knowledge of chemical/physical processes and 

^ the various techniques that are applied to study solar/extra-solar nebular 

i— i chemistry. As such our reference list is limited to works that provide a direct 

discussion of the topic at hand. 

A few aspects of general terminology and background are useful. I will use 
QO riH to refer to the space density (particles per cubic centimeter) of atomic 

f^*) hydrogen, njj 2 for the molecular hydrogen density, and n to the total density 

I"""""" (riH + «h 2 ) ■ Similar terminology will be used for the column density (particles 

per square centimeter) which is denoted by AT. In addition T gas will refer to 
the gas temperature and Td us t for the dust/grain temperature which are not 
always equivalent and have separate effects on the chemistry. At high densities, 
where they are equivalent, I will simply refer to the temperature, T. 

Stars are born in molecular clouds with typical densities of a few thousand 
H2 molecules per cubic centimeter, gas temperatures of ~ 10 — 20 K, and 10 3 — 
10 5 M Q . These clouds exist over scale of tens of parsecs (1 pc = 3 x 10 18 cm), 
%1 but exhibit definite substructure with stars being born in denser (71 > 10 5 

cm~ 3 ) cores with typical sizes of 0.1 pc. The chemistry of these regions is 
dominated by reactions between ions and neutrals, but also by the freeze- 
out of molecules onto the surfaces of cold dust grains. A summary of the 
properties of the early stages can be found in the reviews of Bergin & Tafalla 
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(20071 and di Francesco et al. (20071. Molecular cloud cores are rotating and, 
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upon gravitational collapse, the infalling envelope flattens to a disk. This stage 
is labeled as Class 0. As the forming star begins to eject material the envelope 
is disrupted along the poles and the star/disk system begins the process of 
destroying its natal envelope. In this "Class I" stage the star accretes material 



from the disk and the disk accretes from the envelope (see Adams et al 



1987 Andre et al. 1993 for a discussion of astronomical classification). As 
the envelope dissipates, the system enters the Class II stage wherein material 
accretes onto the star from the disk. The disk surface is directly exposed to 
energetic UV radiation and X-rays from the star, but also the interstellar 
radiation field. Observational systems in the Class II stage are often called T 
Tauri stars. These disks have strong radial and vertical gradients in physical 
properties with most of the mass residing in the middle of the disk, labeled 
as the midplane. It is the midplane that is the site of planet formation. The 
density of the midplane significantly exceeds that of the dense natal core, 
by many orders of magnitude. The midplane is in general, colder than the 
disk upper layers, which can also be called the disk atmosphere. Finally, the 
outer regions of the disk (r > 10 AU) has reduced pressure and temperature, 



but contains most of the mass. Bergin et al. (20071 and Ciesla & Charnley 



(2006) provide recent summaries of disk chemical evolution. This chapter will 
provide greater detail on the methodology and techniques used to explore the 
chemistry of protoplanetary disks and the reader is referred to the previous 
reviews for additional information regarding both observations and theory. 

We will start our study by discussing some key observational results, be- 
ginning with bodies in our solar system and extending to circumstellar disk 
systems at a typical distance of 60-140 pc. We expand our exploration to the 
basic facets of theoretical studies of disk chemistry: thermodynamic equilib- 
rium and chemical kinetics. This is followed by an outline of our knowledge 
of key physical processes. In the final sections, we synthesize these results to 
summarize our theoretical understanding of the chemistry of protoplanetary 
disks 



2 Observational Constraints 

There exists a wealth of data on the chemical evolution of our own solar sys- 
tem regarding the composition of planetary bodies, moons, asteroids, comets, 
and meteorites. In addition, there is a rapidly growing observational sample 
of molecular and atomic transitions in extrasolar protoplanetary disks. The 
current astronomical dataset is limited only to a few objects which have been 
subject to deep searches. However, the upcoming Atacama Large Millimeter 
Array (ALMA) will dramatically expand our current capabilities and offers 
great promise for gains in our understanding of planet formation. This will 
directly complement knowledge gained from the remnants and products of 
planet formation within the solar system. For more detailed views of disk 
chemistry the reader is referred to the recent compilations found in Mete- 
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orites and the Early Solar System II ( Lauretta & McSween 2006 ) , the Trea- 



tise on Geochemistry (Vol. 1, Davis et al. 2005 1, and Protostars and Planets 



V (Reipurth et al. 2007). 



2.1 Planets, Comets, and Meteorites: Tracing the Midplane 

This paper cannot do justice to the decades of study gained in the study of 
bodies within our solar system. Instead we will highlight a few basic facets of 
our understanding and a few examples where the chemical composition is a 
fossil remnant that tracks conditions which existed billions of years ago. 



Planets 



The standard model for the formation of the Earth and other rocky planets is 
that the temperature of the nebula was initially hot enough that all primor- 
dial grains sublimated. As the nebula cooled various species condensed out of 
the hot gas depending on their condensation temperature (defined in § 3.1|. 



The nebula also had strong radial temperature gradients such that different 
species could first condense at larger radii. This model and various issues 



Palme et al. 


(1988 


1 and 


Davis 


(2006 



condensation temperatures of various minerals and ices along with the radii 
where these species potentially condensed in the nebula as estimated by [Lewis] 



(19901. The fraction of rock/ices in planets is also provided on this plot and 
clearly follows this trend. It is a basic astronomical fact that the rocky planets 
reside in the inner nebula, with gas giants in the outer; however, this is sup- 
plemented by Uranus and Neptune having incorporated a greater percentage 
of ices. Thus, basic molecular properties (i.e. temperature of condensation) 
with the nebular thermal structure work in tandem to determine the compo- 
sition of bodies in the nebula to be supplemented by dynamical evolution (e.g., 



Morbidelli et al. 



2000 1 . Other facets of chemical composition of bodies in the 
solar system in relation to nebular thermal structure are discussed by |Lewis" 
(19741, providing some of the basis for the condensation temperature/radii 



placement in Fig. [T| 



Meteorites 



Meteorites are the most directly accessible primitive material in the solar sys- 
tem. Meteorites come in a number of classes, based, in part, on whether the 
material is undifferentiated and therefore more primitive, or part of the mantle 
or core of some larger planetesimal (see Krot et aTj 2005[ ) . The most primitive 
meteorites are the carbonaceous chondrites (labelled as C) . Amongst the vari- 
ous types, the CI carbonaceous chondrite contains a composition that is quite 
similar to the Sun (Palme & Jones 2003[ ). In fact in many instances abun- 



dances from CI chondrites provide a better value for the Solar composition 
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Fig. 1. Top: Plot of temperature of condensation for materials that represent some 
of the major contributors to planetary composition as a function of distance from the 
Sun. This part of the figure is adapted from NASA's Genesis mission educational 
materials. Condensation temperature and estimated position from |Lewis| ( |1990[ ). 
Bottom: Plot of percentage of planetary mass constituted by rocks (black filled 
labels) and ices (white filled labels). Labels represent each of the 8 major planets in 
the solar system. Jupiter, Saturn, Uranus, and Neptune estimates are from |Guillot| 
(1999a I and Guillot (1999b I. Values for the giant planets have larger uncertainties 



than terrestrial planets. 
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Fig. 2. Plot of abundances measured in C I Chondrites against values estimated in 



the Solar Photosphere using compilation of Lodders (20031. 



than spectroscopy of the solar photosphere (Palme & Jones 2003 Lodders 



2003 1 . This is illustrated in Fig. [2] where we plot abundances measured in CI 



chondrites against those estimated in the Sun. 

The exceptions to the general agreement in Fig. [2] are telling. Li is under- 
abundant in the Sun due to convective mixing into layers where Li is destroyed. 
Li is inherited from the interstellar medium and stars begin their evolution 



with abundant Lithium. Primordial Li is gradually destroyed ( Bodenheimer 



1965 Skumanich 19721; thus the presence of Li in a stellar photosphere is 



generally taken as a sign of youth. The most volatile elements have low abun- 
dances in chondrites, this includes the noble gases, along with H, N, C, and to 
some extent O. Unreactive noble gases reside for the most part in the nebular 
gas along with H (which was in molecular form). Nobel gases can be trapped 
inside the lattices of rocky and icy planctcsimals and their abundances pro- 
vide key clues to the evolution of terrestrial planet atmospheres ( |Pepin[ 1991 



Owen fc Bar-Nun[ |1995| |Dauphas| |2003[ ). In the case of C and N, thermo- 



dynamic equilibrium would place these molecules in the form of CO / CH4 or 
N2/NH3 (depending on the pressure/temperature). Thus the temperature in 
the nebula where these rocks formed must have been above the sublimation 
temperature for these molecules. In the case of O 1, the oxygen not consumed 
by refractory material and CO likely resides in the form of gas-phase H 2 0. 
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Fig. 3. Plot of elemental abundances measured in the CM/CO/CV chondrites nor- 
malized to CI shown as a function of the 50% condensation temperature for that 



element. Figure adapted from Alexander et al. (2001 1. Elemental abundances within 



chondrites are taken from Wasson & Kallemeyn ( 1988 1 with condensation temper- 
atures from |Lod dcrs (2003[). 



An additional clue to the chemistry of rock formation in the inner neb- 
ula lies in the relative abundances of elements within the various classes of 
meteorites. Fig. [3] shows a sample of these results plotting the abundances of 
CM, CO, and CV carbonaceous chondrites normalized to CI as a function of 
condensation temperature. For a more expanded plot, the reader is referred to 



Davis ( 2006 1 . This plot demonstrates that elements with higher condensation 
temperatures are more likely to be incorporated into planetesimals. Hence, 
rock formation in the solar nebula was in some sense a volatility controlled 
process. This can involve aspects of both condensation in a cooling nebula 



and potentially evaporation during heating events (Palme et al. 1988 Davis 
20061. 
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Comets 



Because they harbor volatile ices, comets are likely the most pristine rem- 
nants of the formation of the solar system, tracing conditions that existed in 
the outer tens of AU in the solar nebula. Long period comets, with periods 
> 200 yrs, originate in the Oort cloud, while short period comets (< 200 yrs) 
are believed to originate in the Kuiper belt. Oort cloud comets are believed to 
have formed predominantly in the region between Uranus and Neptune and 



scattered to distant orbits via interactions with the giant planets (Fernandez 



1997 Dones et al. 20041. The chemical composition of volatiles in cometary 



coma can be studied via remote sensing techniques with the basic physics sum- 
marized in § |2.2| In general, there are some gross similarities in the abundances 
of volatiles in cometary coma to that of interstellar ices, which are suggestive 



of a direct link (e.g. Bockelee-Morvan et al. 2000 Ehrenfreund et al. 20041. 



In particular water, CO, and CO2, are the most abundant ices in comets and 
in the ISM, with smaller contributions from CH4, NH3, and others. However, 
there is also a wide degree of chemical inhomogeneity in the cometary inven- 



tory. In particular a study of 85 short-period comets by A'Hearn et al. ( 1995 1 



found that a fraction appeared to be depleted of the carbon chain precursors 
that produce C2. Such diversity is not limited to short-period comets as Oort 
cloud comets also exhibit variations in CO, C2H6, and CH4 relative to H2O, 



along with other compositional differences (sec Bockelee-Morvan et al. 2004 



Disanti & Mumma 2008 and references therein). Methanol also stands out 



as another species where variations are found within the Oort cloud popula- 



tion of comets and also short period comets (Disanti & Mumma 20081. These 



differences likely point to systematic changes in the chemistry throughout the 
cometary formation zone perhaps due to a diversity of formation radii for 
cometary nuclei. 

One potential clue to the formation of cometary volatiles lies in the or- 
tho/para ratio of water ice. The spin pairing of the hydrogen atoms for H2O 
leads to two independent forms, para-H^O (spins anti-parallel with total nu- 
clear spin I p = 0) and ortho-H20 (spins parallel with total nuclear spin 
I p — 1). The lowest energy level of para-H20 is ~ 34 K below that of ortho- 
H 2 and if water forms at temperatures below this difference then water will 
preferentially be in the form of para-H20. The ortho/para ratio depends on 
the distribution of population within the energy levels that are characterized 
by three quantum numbers (J, K +1 KJ) used for asymmetric tops (§ 2.2 1. In 



thermodynamic equilibrium the ratio is given by ( |Mumma et"aH|1987| 



o/p 



2I + 1 S{2J + l) e -Eo{J,K + ,K_)/kT 



2I„ 



-E p (J,K + ,K_)/kT ' 



(1) 



1 S(2J+l)e 

and at high temperature is set by the ratio of the spin statistical weights 
(21 + 1), o/p = 3:1. When water forms in the gas via exothermic reactions 
the energy is much greater than the ground state energy difference and the 
ortho/para ratio reflects the 3:1 ratio of statistical weights between these 
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species. When water forms on the surfaces of cold dust grains it is believed 
that the excess energy in the reaction is shared with the grain and the water 
molecules equilibrate to an ortho/para ratio at the grain temperature (e.g. 
|Limbach et al.[ |2006[ ). If the grain temperature is below ~ 50 K then the 
ortho/para ratio will he below 3:1. In Comet P/Halley, Mumma et al. (1987) 
measured an o/p ratio of 2.73 ± 0.17. This is consistent with an equilibrium 
spin temperature of ~ 30 K, which is similar to that measured for water 



and other species (ammonia, methane) in other cometary coma (Kawakita 
et al.| [2006}. | Bonev et ah] ( |2007[ ) provide a summary of water ortho/para 



ratios measured in cometary coma which span a range of spin temperatures 
from ~ 20 — 40 K. One potential interpretation of this ratio is that cometary 
ice formed at temperatures near 30 K, possibly in the interstellar medium, 
although other explanations exist (see Bockelee-Morvan et al. 2004 for a 
discussion) . 



Isotopic Chemistry 



A clear chemical constraint on planet formation is the deuterium enrichment 
observed in the Solar system. In Fig. [4] we provide a compilation of (D/H) 
ratios seen in a variety of Solar System bodies and molecules in the ISM. This 
figure essentially summarizes salient issues as determined by cosmochemical 
and geological studies over the past decades. (1) Gas in the early solar neb- 
ula (proto-Sun, Jupiter) is believed to have a ratio comparable to that seen 
in the main mass reservoir in the natal cloud (i.e., ISM H2). (2) Moving to- 
wards bodies that formed at greater distances, there is a clear trend towards 
deuterium enrichment in trace gas components (Uranus, Neptune). (3) The 
material that formed directly from ices and rocks (comets, IDP, meteorites) 
exhibit large enrichments. (4) The Earth's oceans are enriched in deuterium. 

Based on astronomical studies of the dense (n > 10 5 cm -3 ) cores of molec- 
ular clouds - the sites of star and planet formation - we know that these 
enrichments are readily initiated by a sequence of reactions that are favored 

Indeed 



3.2 



Geiss & Reeves 



1981 



Millar et al. 



1989). 



for T gas < 30 K (§ ^ 
large (D/H) ratios are observed mside mo 

phase and in evaporated ices (Fig. |4| [Ceccarelli et al. 20071. The same cold 



ecular cloud cores, both in the gas- 



(Tgas <j 30 K) non-equilibrium gas chemistry that gives rise to the deuterium 
fractionation in molecular cloud cores is also observed to be active in the outer 
reaches of extrasolar protoplanetary disks ( |Bergin et al.| 2007[ and references 
therein). These inferences hint that the Earth received a contribution to its 
water content from some cold reservoir (e.g. comets, icy asteroids). 

Beyond hydrogen, meteorites also exhibit a number of other isotopic 
anomalies (see Lauretta & McSween 2006 Davis et al.[ 2005[ for a com- 



plete discussion) . Of particular interest are the isotopic enhancement seen in 
carbon, oxygen, and nitrogen which are potentially related to kinetic chemi- 
cal effects early in nebular evolution or perhaps in the cloud that collapsed to 
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Fig. 4. D/H ratios measured in key species in Solar System bodies and molecules 
in the ISM. Plot adapte d from [Messenger et al.| ( |2003[ ) using additional ratios taken 



from the tabulation of Robert ( 2006 1 with references therein. The Jupiter value 



is an in situ measurement from the Galileo probe ( Mahaffy et al. 1998 \ and ISO 



observations (Lellouch et al. 19961, while the protosolar rati o is estimated from 
measurements of J He in the Solar wind ( Geiss & Gloeckler 1998). SMOW - Standard 
Mean Ocean Water - is the reference standard for the Earth's ocean. The ISM value 
is based upon measurements of atomic D/H in clouds within 1 kpc of the Sun ( |Wood| 



et al. 2004 1 
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2.2 Astrophysical Techniques 

Astronomy offers several methods to explore the predominantly molecular 
chemistry of protoplanetary disks. In general, molecules are fantastic probes 
of their environment. In the following we will outline a few basic facets of 
molecular astrophysics and discuss the current status of observations. A listing 
of molecules and transitions that have been detected in protoplanetary disks 
is provided in Table [l] 
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Molecular Astrophysics 

For atomic species the electronic states are discrete, or quantized, and their 
emission spectra has been successfully characterized using the principles of 
quantum mechanics (e.g. Cowley et al. 2000 or an astronomical perspective). 



In comparison to atoms, molecules have more complicated structures where 
rotational and vibrational motions are coupled to electronic states. Similar to 
atoms, the various motions/states are quantized. A more complete description 



of the molecular physics and spectroscopy is provided by Townes fc Schawlow 
(19551 and Herzberg ( |1950[ and subsequent volumes). Evans (19991 provides 
a more descriptive outline regarding the use of molecules as astrophysical 
probes. 



The Born-Oppenheimer approximation (Born & Oppenheimer 19271 al- 



lows the separation of the nuclear and electronic motions in the wave func- 
tion. Thus molecular transitions are characterized by three different energies: 
rotational, vibrational, and electronic. Electronic transitions have energies 
typically ~4 eV (40,000 K) with lines found at visible and ultraviolet (UV) 
wavelengths. Vibrational levels probe energies of ~0.1 eV (1,000 K) and are, 
for the most part, found in the near to mid infrared (~ 2 — 20/im). Rotational 
transitions of the ground vibrational and electronic states have energies below 
0.01 eV (100 K) and emit at millimeter and submillimeter wavelengths. The 
wide range of temperatures and densities present in protoplanetary disks has 
led to detections of electronic, vibrational, and low/high energy rotational 
lines of a variety molecular species (Table [TJ . 

For the bulk of the disk mass, which has temperatures below tens of Kelvin, 
the disk is best traced by molecular rotational lines. Solutions of the wave 
equation for a "rigid" rotating linear molecule yield rotational energy levels: 



1) 



(2) 



where / = /iag is the moment of inertia (/i is the reduced mass, and do 
the bond length, typically ~ 1 A) and J is the angular momentum quantum 
number. The quantity h 2 /2I is redefined as the rotation constant, Bq, which is 
an intrinsic property of a given molecule. For linear molecules the excitation 
is simplified because there is symmetry about all rotational axes and the 
selection rules for electric dipole transitions are AJ ± 1. For these molecules 
the frequency of a rotational transition from level J to level J — 1 is, 



v = Ej{J) - E(J - 1) = 2B J. (3) 

Heavy linear molecules with a large moment of inertia (e.g. HCN, CS, CO) 
have small rotational constants and, in consequence, have more transitions 
at longer millimeter/centimeter wavelengths. Lighter molecules, such as H2 
or HD, have larger energy spacings and emit primarily in the near to far- 
infrared (2 — 115/im). Molecules have 3 degrees of rotational freedom and 
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the simplest case is the linear rigid rotor. Some non-linear molecules have 
a degree of rotational symmetry. For a symmetric top molecule two of the 
principal moments of inertia are equal. In this case the rotational energy 
spacing is charactered by 2 quantum numbers (J,K). Due to their energy 
level spacing and excitation properties, symmetric tops have been used as 



probes of the gas temperature. Prominent examples are NH3 or CH3C2H (Ho 



& Townes 1983 Bergin et al.[ 19941. Non-linear molecules with no symmetry 
are called asymmetric tops. The energy levels for these species are described 
by three quantum numbers (J, K + , K-). Water is the most notable species in 
this category, which also includes most complex organic molecules. 

The spontaneous emission arises from a rotating dipole, which from an 
upper to a lower J state of a linear molecule is given by the Einstein A- 
coefficient: 

S 64 ^ 3 2 ,^ 

Here fid is the permanent electric dipole moment of the molecule. A common 
unit for the dipole moment is the Debye (1 D = 10 -18 esu cm). S is the line 
strength of the transition. For linear rotors S = J u . 

From Eqns.|3]and|4] the methodology of using molecular emission to probe 
the physical environment can be demonstrated. In a disk with density nn 2 , 
molecules are excited through collisions with molecular hydrogen at a typical 
rate hh 2 lui, where the collision rate coefficient, 7„; <~ 1CT 11 cm 3 s^ 1 (see 



Flower 2003. for a thorough discussion of molecular collisions). This leads 
to a critical density required for significant excitation of a given transition, 
which is n cr ~ A u i/~f u i. Using Eqn. |3]and|4]the critical density scales 
Bq/j^J 3 . Of these, the rotation constant and the dipole moment are molecular 
properties, while J depends on the observed transition. Heavy molecules with 
large dipole moments, due to charge disparities, have higher critical densities 
and are excited in denser gas. Thus HCN (fid — 2.98 D) and H2O (fid = 
1.85 D) are tracers of dense (rt H2 > 10 5 cm -3 ) gas, and CO (fid = 0.11 D) 
traces lower density material. Higher J transitions probe successively denser 
(and warmer layers), presumably at the disk surface, or radially closer to the 
star. 

Molecules with identical nuclei have no permanent electric dipole, but in 
some instances can emit via weaker quadrupole (H 2 ) or magnetic dipole (O2) 
transitions. In the case of H2, which is the dominant gas constituent, the first 
rotational level ( J = 2 — for a quadrupole) is 510 K above the ground state. 
The population in this state is cx exp(— E u / 'kT gas ) , which hinders H2 emission 
from the cold (T gas ^10 — 30 K) midplane. Hence, similar to the ISM, other 
molecules with dipole moments are used as proxies to trace the molecular gas. 
These species arc generally heavier molecules comprised of H, C, O, and N 
(typically CO). 

Vibrational modes are also commonly detected in warm and dense disk 
systems. Molecules with N atoms have 3N degrees of freedom. In a linear 
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molecule 2 degrees are rotational and 3 translational. Thus there are 3N — 
5 degrees of freedom (3N — 6 for non-linear). In the case of H 2 there are 
3 modes, two stretching and one bending. Vibrational modes are quantized. 
However, the stretching frequency can be roughly approximated in a model 
where the two atoms are attached via a spring. Thus using Hooke's law, where 
the vibration of spring depends on the reduced mass Qit r ) and force constant 
k (the bond), 




The force constant for bonds is ~ 5 x 10 5 Dynes cm -1 for single bonds, 10 
xlO 5 Dynes cm -1 for double bonds, 15 xlO 5 Dynes cm -1 for triple bonds. 
The energy of vibration depends on the vibrational quantum number, n: 

Kib = (n+^jhu . (6) 

Thus heavier molecules have vibrational modes at shorter wavelengths 
when compared to lighter molecules. Each vibrational state has a spectrum 
of rotational states called ro-vibrational transitions. The frequency of these 
transitions from an upper state ' to a lower state " is given by (ignoring some 
higher order terms): 

v = v Q + B V ,J'(J' + 1) - B'1„J{J" + 1). (7) 

B v contains the rotational constant along with a correction due to the inter- 
action between vibration and rotation (i.e. the intermolecular distance can 
change). The selection rule is A J = ±1 so there are 2 cases J' = J" + 1, 
the R branch, and J' = J" — 1, the P branch. These are labelled as R(J) or 
P( J) where J is the lower state. In some instances AJ — is allowed, which 
is called the Q branch. Quadrupole transitions Z\J±2 are labeled as S(J) for 
AJ = 2 and 0( J) for AJ = -2. 



Brief Summary of Observational Results 



Some features of the observational picture become clear in Table [T] Molecular 
transitions that probe higher temperature gas have been found to predomi- 
nantly probe the innermost regions of the disk. This is inferred via spatially 
unresolved but spectrally resolved line profiles, assuming that the line broad- 
ening is due to Keplerian rotation. Spatially resolved observations of the lower 
energy rotational transitions trace the outer disk and confirm that protoplan- 



etary disks are in Keplerian rotation ( jKoerner fc Sargent 1995 Simon et al. 
2000[ ) and that the temperature structure has sharp radial gradients (as ex- 



pected). In addition, analysis of high- J molecular transitions by van Zadelhoff 



[ 
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et al. ( 2001 1 and Aikawa et al. ( 2002 1 confirms that the surface of the irradi- 
ated disk is warmer than the midplane, as suggested by analysis of the dust 
spectral energy distribution ( Calvet et al.|[l99T| Chiang & Goldreich 19971. 
It is important to state that current sub/milli-meter wave observational facil- 
ities are limited by resolution and sensitivity and, therefore, only probe size 
scales of tens of AU at the nearest star-forming regions (e.g. TW Hya at 65 
pc, and Taurus at 140 pc). This limitation will be significantly reduced with 
the advent of the ALMA array (and JWST) in the coming decade and we can 
expect the listing in Table [T] to drastically increase. 



Table 1. A Sample of Current Astrophysical Probes 



Species 


A(/im) 


Transition 


E u (K) 




Radius Probed 


Notcs a 


H 2 


0.10 


- 0.15 


Lyman- Werner bands 


10° 


r 


< 1 AU 


(1) 


H 2 




2.12 


V = 1 - S(0) 


6471 


r 


~ 10 - 40 AU 


(2) 


CO 




2.23 


v = 2 - 


6300 


r 


~ 0.05 - 0.3 AU 


(4) 


H 2 




~ 2.9 


V 3 = 1 - 


5000 - 10000 


r 


~ 1 AU 


(7) 


OH 




~ 3 


V = 1 — P branch 


> 5000 


r 


~ 1 AU 


(7) 


CO 




4.6 


V = 1-0 


3000 


r 


~< 0.1 - 2 AU 


(5) 


H 2 


8.0 - 


17.0 


v = 0-0 S(l), S(2), S(4) 


1015 - 3474 


r 


~ 10 - 40 AU 


(3) 


H 2 


10 


- 30 


J > 4 


> 500 


r 


~ 1 - 2 AU 


(6) 


C 2 H 2 




' 13.7 


i?5 - 1 - Q branch 


1000 


7* 


~ 1 AU 


(8) 


HCN 




~ 14 


D 2 — 1 — Q branch 


1000 


r 


~ 1 AU 


(8) 


co 2 




14.98 


V2 — 1 — Q branch 


1000 


r 


~ 1 AU 


(8) 


No II 




12.81 


P3/2 — Pl/2 


1100 


r 


~ 0.1 AU b 


(9) 


CO 


460 - 


2600 


6 - 5, 3 - 2, 2 - 1,1-0 


5 - 116 


r 


> 20 AU° 


(10) 


HCO+ 


1000 - 


3300 


3 - 2. 1 - 


5-25 


r 


> 20 AU° 


(11) 


cs 


1000 - 


3000 


2 - 1,3 - 2, 5-4 


5-30 


r 


> 20 AU C 


(12) 


N 2 H+ 




3220 


1-0 


5 


r 


> 20 AU° 


(13) 


H 2 CO 


1400 - 


2000 


3i3 — 2i 2 ,2i 2 — ln,3i 2 — 2n 


20 - 32 


r 


> 20 AU C 


(14) 


CN 


1000 - 


2500 


3-2.2-1 


5-30 


r 


> 20 AU° 


(15) 


HCN 


850 - 


3300 


4 - 3, 2 - 1,1-0 


5-40 


r 


> 20 AU C 


(16) 


HNC 




3300 


1 - 


5 


r 


> 20 AU C 


(17) 


H 2 D + 




805 


lio — In 


104 


r 


> 20 AU C 


(18) 


DCO+ 


830 - 


1400 


5-4,3-2 


20 - 50 


r 


> 20 AU C 


(19) 


DCN 




1381 


3-2 


20 


r 


> 20 AU C 


(20) 



Refer ences : (l)|Herczeg et al!] l |2002), (2) |Bary et al.| <2003| |2008), (3) dBitner 



2007>, (4) lUarr et al.| 119931, (5) I Najita et al.H2UU3| ; |nrittain et al.| jn007[, ( t> ) 
fe IN ajitalJaUUShlBalyk et al.H2008t, (7)|HaIyk et al.| (SSSbf, (« ) |Lahuis et all <2DD 
fe Najital 12008i, (9) IKspaillat et al.H2007|; ILahuis et al.H2007[; IHerczeg et 



C^ar r_&: INaj 

71, (10)|Uutrey"et al 



( 1996 ) 



astner et al.| U997II; |Qi| H200l[ 
~ n^7l ;l^astner et al.| 



et a 

van Zadclhoff ct a 
19971; Ivan Zadelhott 



2UUH;|Qi et al.H2006l, (ll ) |Uutrey et a 
|e t al.|j200H;|gi et al.||2UU«>, (12)|L)utrey et al.H1997>, (13)|Uutrey et al.HTg9"7||200' 

14) |Dutrey et al.| 11997>; |Thi et al.H2UU4|, ( 15) IDutrey et al.| jl997|l; jkastner et al. 
l|1997t;|van Zadelhotf et al.|1200H, (16 ) |Dutrey et al.| U997II; IKastner et al.| l|1997(; |va n| 
k adelkott et al.|12UUH;|(ji et al.|H2008t, (17)|Uutrey et al.| 119971; IKastner et al. |l |199 
a) |(Jeccarelli et aTf]2004( |van Uishoeck et al.|q2003[ ; |Qi et al.| 12008^ 7(2071 

|et al'H^O^ 

b If the [Ne II] emission arises from a photoev aporative wind then the emission can arise 
from greater distances jHerczeg et al.| |2007fl . 

c It is important to note that many ot these species will have rotational emission inside 
20 AU, particularly in the high-J transitions. However, the observations arc currently 
limited by the spatial resolution, which will be overcome to a large extent by ALMA. 



The chemistry of protoplanetary disks is slowly being cataloged. Molecu- 
lar emission depends on physical parameters (density, temperature, velocity 
field) and the molecular abundance. The radial distribution of density and 
temperature is estimated via thermal dust emission (e.g. Beckwith et al. 
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Fig. 5. Infrared spectrum of the classical T Tauri star AA Tau taken by the Spitzer 
Space Telescope. Vibrational modes of C2H2, HCN, and CO2 along with a transi- 
tion of [Ne II] are labelled. Diamonds mark detections OH rotational transitions. 
Numerous rotational transitions of water vapor (not marked) are spread throughout 
the spectrum. Figure and discussion published by Carr & Najita ( 2008[ ). 



19901. Power-law fits provide £(r) oc r~ p and T(r) oc r~ q with p = — 1 
and q = 0.5 — 0.75. Temperatures are estimated to be ~ 100 — 200 K at the 
midplane at 1 AU with surface densities of ~ 500 — 1500 g cm -2 . Modern 
models now take into account the disk vertical structure (§ 4.1 Calvet, this 



volume; Calvet et al. 1991 Chiang fc Goldreich 1997 D'Alessio et al. 19981 



Thus the determination of chemical composition requires untangling the phys- 
ical and chemical structure, both of which can vary in the radial and vertical 
direction. Abundances were first calculated using local thermodynamic equi- 
librium in a disk with radial variations in density and temperature, but now 
adopt more sophisticated models capable of exploring 2 dimensional parame- 



ter space (Monte-Carlo, Accelerated Lambda Iteration: Bernes 1979 Rybicki 
& Hummer 



1991) 



With these techniques, a simple understanding of the chemistry is emerg- 
ing. One underlying result is that the molecular abundances are generally re- 
duced when compared to those measured in the interstellar medium ( |Dutrey| 
et al.||1997||Bergin et al] |2007| [Dutrey et al.[|2007a[). This is due to the freeze- 
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out of molecules in the dense cold midplane, which will be discussed later in 
this Chapter. Closer to the star there exists species specific "snow-lines" where 
molecular ices evaporate. In this regard, the presence of hot abundant water in 



the inner (< 2 AU) disk has been confirmed in some systems (Carr & Najita 
2008 ) . A sample of these results is shown in Fig. |5j where a 



2008 Salyk et al. 



Spitzer spectrum of a T Tauri star is shown. This spectrum reveals not only 
the presence of water vapor inside the snow-line, but also abundant pre-biotic 
organics (see also Lahuis et al. 2006). It is also becoming possible to resolve 
chemical structure within protoplanetary systems (e.g. Dutrey et al. 2007a I. 
A sample of what is currently possible is shown in Fig. [6] where we show the 
integrated emission distributions of HCO + , DCO + , HCN, and DCN. As can 
be seen the emission of deuterium bearing DCO + is off-center when compared 
to HCO + , which is indicative of an active deuterium chemistry, at least for 
that species. Finally, in the case that most species containing heavy elements 
(with dipole moments) are frozen onto grains in the dense midplane then two 



species, H2D" 1 " and D2H+, will remain viable probes (Ceccarelli et al. 20041 



Some complex organics have been detected in the mm-wave and a summary 
of disk organic chemistry is given by Henning fc Semenov| ( |2008 1 . 

One key factor for modeling disk chemistry is the dissipation timcscale of 
the gas rich disk. The molecular detections listed in Table[T]are predominantly 
from systems with ages of ~ 1 — 3 Myr, demonstrating that the gas rich disk 
clearly exists for at least a few Myr. In addition, some gas-rich systems, such 
as TW Hya, are observed at 10 Myr. There is growing evidence of increased 
dust evolution (i.e. grain growth and possible planet formation) in the inner 
disk (r < 10 — 30 AU) of some systems (Calvet, this volume; Najita et al 



2007 1 . Systematic surveys of disk dust emission in stellar clusters has provided 



some of the strongest limits on the evolutionary timescales for solids. Based 
on these studies it is estimated that the timescale for grain growth, settling, 
and potentially planctcsimal formation, is ~ 1 — 10 Myr (e.g. Haisch et al. 



2001 Cieza 2008 1. The disk gas accretion rate decays over similar timescales, 



although there is a large scatter in the data ( Muzerolle et al.[ 2000). Beyond 
disk accretion, direct observations of gas dissipation is limited by our observa- 
tional capabilities. However, a Spitzer search for gas emission lines in systems 
spanning a wide range of ages has suggested that the gas disk dissipates on 
similar timescales as seen for dust evolution (i.e. ~ 1 — 10 Myr; Pascucci et al 



2006 ) . It is worth noting that the available data, while showing an general de- 



cline with age, exhibits a large degree of scatter. Thus there are systems that 
have lost their disk in ~ 1 Myr, while others exhibit rich gaseous disks at 



10 Myr (Cieza et al. 2007) 



3 Equilibrium and Kinetic Chemistry 



In the following we will outline the basic methods and chemical processes that 
are generally used to examine chemistry in protoplanetary disks. In the past, 
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Fig. 6. Figures showing sample spectra and integrated emission maps of the J = 3 — 
2 transition of HCO+ , DCO+ , HCN, and DCN . Observations from the Submillimeter 
Array and figure published by Qi et al. (2008). 
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the primary method to explore disk chemistry involved assuming a gas-solid 
mixture in fixed pressure and temperature with solar composition and de- 
termining resulting composition in thermodynamic equilibrium. This method 
has been quite successful in explaining various trends observed in meteorites 



and planetary bodies. For a summary of these efforts see Prinn ( 1993 1; Davis 



& Richter ( 2004 and references therein) . More recently, in part due to obser- 
vations of an active chemistry on the disk surface, it has been recognized that 
for much of the disk mass the temperatures and pressures are too low for the 
gas to reach equilibrium and the gas kinetic chemistry must be considered. 
In the following we will describe the basic techniques for both of these cases, 
with a brief discussion regarding the regimes where each is relevant. 

3.1 Thermodynamic Equilibrium 

Thermodynamic equilibrium is defined as the state of maximum entropy or, 
the state of minimum Gibb's free energy of the system. It is the state towards 
which all chemical systems are being driven if one waited "forever and a day" . 
A key aspect of these calculations is that the final composition depends on the 
pressure, temperature, and elemental composition, but not the initial chemical 
state. For the following the author is grateful to C. Cowley for numerous 



discussions and providing an outline of the methodology (see Cowley 1995 1 . 

The simplest method involves using the equilibrium constants. For every 
reaction an equilibrium can be defined. For example in the following reaction: 

2H + O # H 2 (8) 

the ratio of partial pressures (px = nxkT) can be determined from the equi- 
librium constant itp(H 2 0) which is a tabulated quantity (c.f. JANAF tables 
from NIST). The equilibrium constant can be expressed as: 



K P {W) = g^f^ = ex P [-AG°(R 2 0)/RT}. (9) 



p(H 2 Q) 
P 2 (H) P (0) 

Where R is the gas constant and the AG is the Gibbs free energy, which for 
this reaction is, 



Z\G°(H 2 0) = Z\G}(H 2 0) - 2Z\G}(H) - AG° f (0). (10) 

The o superscript denotes that the substance is in its standard state given the 
temperature (i.e. liquid, gas, solid) and subscript / refers to "of formation" . 
In a system with multiple species at fixed pressure and temperature we can 
use the equilibrium constants to determine the composition. For an example, 
lets take the following system: 

CO^C + O (11) 
CN^C + N. (12) 
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We know the following: 



K C o(T) = (13) 
Pco 

KMT) = ™. (14) 

PCN 

We can now define "fictitious pressures", Pq, Po, an d Pn which can be ex- 
pressed, 



Po=Po+ Pco =Po + -r? — (15) 
d , , P NPc 

ftr = Pn + Pcn =Pn+ -rz — (16) 

D i i , PcPn . PcPo 

Pc =Pc+Pcn+Pco =Pc + ^ 1- — (17) 

and the pressure (mass) balance equation, 

P = Pc + Pn + Po + Pco + Pcn (18) 

From these you can solve for the partial pressures and composition. 

There are some other more commonly used methods in the literature, 
which we will outline below. The Gibb's free energy (G) is defined by: 

G = H-TS = E + PV-TS (19) 

where H is the enthalpy, E the energy, S the entropy. P, T, and V, are, 
respectively, the pressure, temperature, and volume. From this the change in 
energy is, 

dG = dE + PdV + VdP - TdS - SdT, (20) 
which, using the 1st law of thermodynamics, reduces to, 

dG=VdP- SdT. (21) 
If this is an ideal gas than PV = nRT, and 

(dG) T - VdP = ^^dP, (22) 

G B -G A = J"" ^j-dP = nRTln (^j (23) 
In the case of departures from an ideal gas, then: 



The Chemical Evolution of Protoplanetary Disks 



19 



Gb — G a — 



PA 



V- 1 ^^dP + nRTln(^j=nRTln(^\ (24) 



where / is the fugacity, which is essentially a pressure term for a non-ideal 
(real) gas. It is equivalent to the partial pressure if the gas is ideal. It is now 
common to re-write this in terms of some known or measured state (e.g. 1 
atm where / = /o): 



G B = G° 



iRTh 



= G° + nRTln{a) 



(25) 



where a is the activity which is defined as f/fo- If we have an ensemble of 
atoms and molecules we can write out the total Gibbs free energy: 



G = £rii. p (AGi.p + RTlna,i p ) 



(26) 



where n is the number of moles of species i in phase p (either pure condensed, 
condensed, or gaseous). The final equation is for mass balance normalized to 
the assumed initial composition. Given these sets of equation the minimized 
value of the Gibb's free energy determines the equilibrium composition. Ther- 



modynamic data are available via the JANAF tables (Chase 1998), although 



the tables are not complete. Some useful references that more outline this 
procedure and more directly discuss issues regarding thermodynamic data 
are IBurrows & Sharpl (119991) andlLodders & Fegleyl (12002 1. 



Condensation 



The condensation temperature of an element as defined by Richardson & 
McSween (1989) is "the temperature at which the most refractory solid phase 



containing that clement becomes stable relative to a gas of solar composition" . 
A full discussion of this topic is beyond the scope of this chapter and the 



reader is referred to Grossman ( 1972 1, Richardson & McSween ( 1989 ), Lodders 
(2003J) (and references therein). However, as discussed in §|2j these calculations 
lie close to the heart of understanding the composition of some meteorites and 
as such warrant a brief discussion. The essential idea is that as the nebula gas 
cools (at constant pressure) refractory minerals condense from the gas and this 
sequence depends on the condensation temperature of the various potential 
refractory minerals. When a given element condenses one must pay attention 
to mass balance as some (or all) of its abundance is now unavailable for the gas 
in equilibrium. As an example, Corundum (AI2O3) is one the first minerals 
to condense. This is consistent with the detection of this minerals in Ca-Al 
rich inclusions in chondrites and the fact that Ca-Al inclusions provide the 
oldest age measured in the Solar System ( Amelin et al.| [2002 ) . An excellent 
further discussion of this topic is found in the references above and also in 



Davis (20061. 
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Fischer- Tropsch Catalysis 



Fischer- Tropsch catalysis is the name for a subset of high temperature re- 
actions wherein CO (and CO2) with H2 is converted to methane and other 
hydrocarbons on the surfaces of transition metals, such as Fe and Ni. If Fe/Ni 
grains are present then Fischer- Tropsch catalysis is an effective method to cre- 
ate hydrocarbons from gas-phase CO at temperatures where the reaction is 
inhibited due to high activation energies (T < 2000 K). The conversion of CO 
to organics via this process has been explored by numerous authors and the 



reader is referred to Fegley & Prinn ( 1989 1, Kress & Tielcns (2001 ) and Sekine 



et al. (2005 1, and references therein. The efficiency of Fischer- Tropsch catalysis 
is highly dependent on the gas temperature and pressure. The efficiency mea- 
sured at pressures relevant to giant planet subnebulsQ (0.09-0.53 bar) peaks 
at ~550 K with reduced efficiency at either lower and higher temperatures 



( Sekine et al. 2005 1 , but see also Llorca fe Casanova| ([2000 ) for measurements 
at lower pressures. At higher temperatures the surface of the catalyst is poi- 
soned by the conversion of surface carbide to unreactive graphite. At lower 
temperatures the efficiency will ultimately be limited by the incorporation of 



Fe into silicates, the formation of Fes04 coatings (Prinn & Fegley 19891, or 



coatings by organic and water ice (jSekine et al. 2005 1 . It is also worth not 



ing that ammonia could also be formed from N2 via a similar process called 
Habcr-Bosch catalysis. 



Clathrate Hydrates 

Lunine & Stevenson ( 19851) provided a detailed analysis of the thermody- 



namics and kinetics of clathrate hydrates at relevant conditions for the solar 
nebula. A clathrate hydrate is a water ice crystal in which the lattice is com- 
posed of water molecules with open cavities that are stabilized by the inclu- 
sion of small molecules of other chemical species. The inclusion of additional 



molecules in the lattice is actually key to the stability of the clathrate ( Buf- 



fett 


2000 


Devlin 


2001) 



open then cither the regular structured lattice of crystalline ice or irregular 



structured lattice of amorphous ice (see, e.g. Tanaka & Koga 20061. 

The relevance of these compounds for nebula chemistry is that they provide 
a means of trapping volatile species (noble gases, CO, N 2 ) into ices at tem- 
peratures well above their nominal sublimation temperature. Thus volatiles 
can be provided to forming planets to be later observed in their atmospheres. 



Lunine & Stevenson (19851 suggest that clathrate hydrates can be important, 



but only if fresh ice is exposed to the surrounding gas via frequent collisions or 



perhaps if ices evaporate and condense upon an accretion shock (e.g., Lunine 



et al. 1991) 



1 Giant planet subnebula refers to the collapsing pocket of the surrounding disk 
gas out of which the giant planet forms. Both the gas/ice giant and its satellites 
form within the subnebula. 
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Kinetic Inhibition 



It was recognized that the conditions in the solar nebula may not always be in 
thermodynamic equilibrium. In particular, the timescale of nebular evolution 
(e.g. the gas lifetime) could be shorter than reaction timescales and therefore 
equilibrium would not be attained. Alternately, the dynamical mixing of gases 
from colder regions within the nebula to the warmer central regions could pre- 
clude the formation of certain species which would otherwise be favored given 
local conditions. Thus knowledge of the reaction kinetics is required. Cases 
where kinetics limit the outcome of thermodynamic equilibrium are labeled as 



kinetic inhibition by Lewis & Prinn (1980). It is also referred to as quenching 
in the literature. As an example, in equilibrium CH4 should dominate over 
CO in the inner solar nebula and also within giant planet subnebula. How- 
ever, the chemical timescale for CO conversion is ~ 10 9 s at T gas ~ 1000 K 
and longer at lower temperatures (Lewis & Prinn 19801. Estimates suggest 



that dynamical timescales are shorter, and the zone where CH4 dominates 
therefore shrinks considerably (e.g. Fegley & Prinn[ |1989[ IMousis & Alibert 



2006[ ). In general, these models derive the equilibrium results and then explore 
whether kinetics would change the predicted results (see, Smith 1998 for a 
discussion) . 



3.2 Chemical Kinetics 



For much of the disk mass the density and temperature are low enough that 
the gas does not have sufficient time to reach equilibrium and one needs to 
explore the time-dependence of the chemistry. The exact temperature and 
density (and various combinations) where equilibrium calculations are not 
valid is presently unclear, and we include a brief discussion of this in § 3.3 



However, it is certain that within regions of the disk where the gas temperature 
is below 100 K (which is most of the disk mass) that the chemical kinetics 
should be considered. 

Fortunately, observations (see § [2]) appear to suggest that the observed 
chemistry on tens of AU scales is similar to that seen in the dense [n > 10 6 
cm" 3 ) regions of the interstellar medium exposed to energetic radiation fields. 



Models developed for application in this regime (Hcrbst 



( 


Herbst 


1995|) have been 


Aikawa et al. |1997 


Bauer 



that the cold (T gas < 100 K) temperatures require exothermic reactions, while 
the low (n ~ 10 3 — 10 6 cm" 3 ) density limits reactions to two bodies. At these 
energies typical reactions between two neutrals are endothermic or have an 
activation barrier. Therefore, exothermic reactions between ions and neutrals 



are believed to dominate (Herbst & Klemperer 1973 Watson 1976 1. However, 



given the strong radial gradients in density and temperature, it is certain that 
in some instances reactions with barriers are important and 3-body reactions 



can also be activated (Willacy et al. 1998 Tscharnuter & Gail 20071. 
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In the disk context the temperature is often low enough that water and 
other more volatile species will freeze (or even be created in situ) on grain sur- 
faces. Thus the inclusion gas-grain interaction is a key piece of any successful 
disk chemistry model. Below we will first describe the basic gas-phase chem- 
istry and how networks are created and solved. Then we extend the discussion 
to include grain chemistry. 



Gas-Phase Chemistry 

In Table [2] we provide a listing of common chemical reaction types and typical 
rates. To be active the gas-phase chemistry requires a source of ionization, 
w hich in the disk can cither be cosmic rays, X-rays, or active radionuclides 
(§ 4.2 1. The ionization of H2 ultimately produces H^", which powers the ion- 
molecule chemistry. A general feature of this chemistry is the subsequent 
creation of complex molecular ions which undergo dissociative recombination 
to create both simple and complex molecules. To illustrate the various features 
of the chemistry we will examine the chemistry of water which has significance 
in the disk context. 



Table 2. Typical Reaction Types and Coefficients 



Reaction Type 



Example 



Typical Rate Coefficient 



Ion-Neutral 


H+ + CO- 


-» HCO++ H 2 r 


- 10" 


" 9 cm 3 s" 1 


Radiative Association a 


C++ H 2 -» CH+ + hv r 


- 10" 


-17 


cm 3 s 1 


Neutral Exchange 


O + OH - 


O a + H 


- 10" 


-11 


- 10~ 10 cm 3 s 


Charge Transfer 


C++ s -> 


S++ c 


- 10" 


" 9 cm 3 s" 1 


Radiative Recombination 3 


C+ + e~ - 


-» C +hu 


- 10" 


-12 


cm 3 s 1 


Dissociative Recombination 3 


H 3 0+ +e 


- -> H 2 + H r 


- 10" 


" 6 cm 3 s" 1 


Photoionization 


C +hv -f 


C+ + e" 


- 10" 


" 9 e 


~" /j ^ v cm 3 s- 1 


Photodissociation 


H 2 Q +hv 


-c O + OH r 


- 10" 


" 9 e 


cm 3 s" 1 



Sample rate coefficients derived at a gas temperature of 40 K using the 
UMIST database fWoodall et al.||2007a|). 



b The depth dependence of photoionization rates and photodissociation 
rates is often parameterized by e ,Av . In this form 7 includes factors 
such as the shape and strength of the radiation field, the absorption 
and scattering due to dust, and the wavelength dependence of the pho- 
toionization or photodissociation cross-sections. For greater discussion please 



see van Dishoeck| ( |1988[ ), |Roberge et al. ([1991]) , and van Dishoeck et al. (2006) 
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The oxygen chemistry begins with the ionization of H2 and the formation 
of the trihydrogen ion (H3"). will react with O and a sequence of rapid re- 
actions produces H30 + . The next step involves the dissociative recombination 
of H30 + with electrons: 



H 3 CT 



H + H 2 
OH + H 2 
OH + 2H 
O + H + H 2 



fx 
f 2 
fs 
f 4 



(27) 



where fi-4 are the branching ratios^] Via this sequence of reactions water 
vapor is then created. The primary gas-phase destruction pathway in shielded 
gas is via a reaction with He + (another ionization product); at the disk surface 
photodissociation also plays a key role. It should be stated that the example 
above provides only the primary formation and destruction pathways under 
typical conditions (i.e. where ionizing agents are available). In addition, water 
participates in numerous other reactions both as reactant and product and 
some water can be created via other pathways, albeit at reduced levels (e.g. 



Bergin et al. 2000 1 



At temperatures below a few hundred Kelvins ion-neutral chemistry dom- 
inates. Above this temperature, two neutral-neutral reactions rapidly trans- 



form all elemental oxygen in the molecular gas into water vapor (Elitzur & de 



Jong 1978 Wagner & Graff 19871 



+ H 2 -> OH + H (E a = 3160 K), (28) 
OH + H2 -> H 2 + H (E a = 1660 K). (29) 

This mechanism (along with evaporation of water ice) will keep water vapor 
the dominant oxygen component in the gas within a few AU of the forming 
star. 

To model these reactions one creates a series of ordinary differential equa- 
tions, one for each species included, accounting for both formation rates and 
destruction rates: dn/dt — j ' or motion — destruction. For water (excluding the 
numerous reactions not mentioned above) the equation is as follows: 



dn(R 2 0) 
dt 



n(H 3 0+)n(e-)a r + n(OH)n(H 2 )k 1 



l (H 2 O)[n(He + )k2 + k 7 ,A v=0 e 



-X.7A V 



(30) 



Two measurements using a storage ring give roughly consistent results ( |Jensen| 



et al. 2000 Vejby-Christensen et al. 



ment: f x = 0.25, f 2 = 0.14, f 3 = 0.60, and / 4 



(flowing afterglow), Williams et al. (19961 find /1 
f 4 = 0.30. 



19971. We provide here the latest measure- 
0.01. Using a different technique 
0.05, h = 0.36, f 3 = 0.29, 
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In the equation above a r is the recombination rate of HaO" 1 " (Jensen et al 
2000 1 and fc 7 ,A„=o is the unshielded photodissociation rate. Values for this rate 



and others can be found in either of the two available astrochemical databases 
(one originally maintained by E. Herbst at Ohio State or the other at UMIST: 



Woodall et al. 2007b I. These rate networks have thousands of reactions oc- 



curring over a wide-range of timescales and thus it is a stiff system of ordinary 
differential equation (ODE). A number of ODE solvers have been developed 
to solve these systems with the output of composition as a function of time. 
A summary of some of the available codes for astrochemical applications is 



given by Nejad (20051. 



Gas-Grain Chemistry 

In the dense (n > 10 5 cm" 3 ), cold (T < 20 K) centers of condensed molecular 
cores, which are the stellar birth sides, there is ample evidence for the freeze- 
out of many neutral species onto grain surfaces ( |Bergin fc Tafalla] |2007[ ) . This 
even includes the volatile CO molecule that only freezes onto grains when 
temperatures are below 20 K (jCollings et al. 20031, as opposed to water 

The 



2001) 



which freezes at 110 K for interstellar pressures (Fraser et al 
outer disk midplane (r > 20 — 40 AU) is believed to be similarly cold with 
densities that are several orders of magnitude higher than in pre-stellar cores 
where volatile CO freeze-out is readily detected. Thus to model disk kinetic 
chemistry it is necessary to model both the deposition of species onto grain 
surfaces in tandem the mechanisms which desorb frozen ices. A number of 
such desorption mechanisms have been identified in the literature including 
thermal evaporation, X-ray desorption, and UV photo-desorption. Each of 
these will be described below and are generally included in the rate equations 



(Eqn. 30 1 as source or loss terms with rate coefficients given below. 



Gas-phase Freeze-out: Atoms and molecules freeze onto grain surfaces or ad- 



sorb with a rate of k/ 



n g a g vS s 



Here a g is the grain cross-section, 



is the mean velocity of the Maxwcllian distribution of gaseous particles, n g 
the space density of grains, and S is the sticking coefficient (i.e. how often a 
species will remain on the grain upon impact). Sticking coefficients for molec- 
ular species are generally thought to be unity at low T which is supported 



by theoretic al ( Burke fc Hollenbach] 19831 and laboratory work (Bisschop 
et al. 2006). For an interstellar grain size distrib utiorj^] with a lowe r cutoff of 



2 x 10 21 n cm 1 (see discussion in Irlolienbach et al. 20081. In 



20 A, n g a g 

disk systems the process of grain coagulation and settling will locally reduce 
the grain surface area and this value is only a starting point for potential 
situations. Assuming this surface area, the freeze-out timescale, tj , is: 



In the interstellar medium grains are inferred to exist with a size distribution 



that follows a power-law of dn(a)/da oc a , with a as the grain radius (Mathis 



et al. 



1977). 
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a) Depletion Timescale b) CO ice Evaporation Timescale 




Fig. 7. (a) CO freeze-out timescale as a function of density of molecular hydrogen, 
(b) Thermal evaporation timescale for CO molecules as a function of dust temper- 
ature. 



r ,«3 /10 5 cm- 3 \ {20K\ ' r- 
, /o = 5xlOVs(— — J VA (31) 

where A is the molecule mass in atomic mass units. Fig. [^demonstrates that 
within the disk the freeze-out time for CO is on the order of tens of years 
or below. Thus, without effective means to remove molecules from the grain 
surface, most heavy molecules freeze onto grain surfaces. In the astronomical 
literature this timescale is often referred to as the depletion timescale. Because 
molecular abundances can "deplete" by freeze-out or gas-phase destruction, 
we adopt here the more descriptive terminology. 

Grain Binding Energies: 

The removal of ices from the refractory grain core requires overcoming 
the strength of the bond to the surface. Because low temperatures preclude 
the breaking of chemical bonds, molecules are not expected to be chemically 
bound to grain surfaces (although chemical bonds may be important at higher 



temperatures; Cazaux fc Tielensj 2004 1. Instead, the approaching molecule has 



an induced dipole interaction with the grain surface or mantle and are bound 
through weaker Van der Waals-London interactions called physical adsorp- 
tion ( |Kittel[ [1976 ). Van der Waals interactions are proportional to the prod- 



uct of the polarizabilities (a) of the molecule and the nearby surface species 
{Eb oc a mo ia sur f ace ). Because of this property binding energies on a silicate 
surface have been calculated using the extensive literature on the physical ad- 
sorption of species onto graphite. Accounting for differences in polarizability 



these can be scaled to silicates (Allen & Robinson 1977). Hasegawa et al. 
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Table 3. Binding Energies of Important Ices a 

Species Eb/k Reference 
H 2 5800 Frascr ct all 2001) 



CO 


850 


Pollings et 1I^1[2TO3| 


N 2 


800 


Bisschop etlu^(|2TIu6| 


NH 3 


2800 


Brown fc B^Ima||2TO7 



CH3OH 5000 iBrown k Bolina (2007) 

1 — 

Values are from laboratory work only 

using pure ices. 



( 1992 1 re-examined this issue and provide a sample list of estimated values 
that have been widely used. These estimates need to be revised for cases 
where actual experimental information exists. Measurements have been now 
been performed for a number of key species, with a listing given in Table [3] 



Thermal Evaporation: 

The rate of evaporation from a grain with a temperature of Tdust , is given 
by the Polanyi-Wigner relation ( Tielens[ 20051: 



kevap.i = VQ,iexp(-Ei/kT dust ) (s x ). 
vq is the vibrational frequency of molecule i in the potential well. 



(32) 



1.6 x 10 1 V(£iA)("W™i) 



(33) 



with Ei, i the binding energy of species i. Following Hollcnbach et al. (2008) 



the sublimation or freeze-out temperature can be derived by setting the flux 
of desorbing molecules (F t d,i) equal to the flux of species adsorbing from the 
gas. Thus, 



Fi 



td.i 



N k 



(34) 
~ 2 ), rii 



where N s i is the number of adsorption sites per cm 2 (N s i ~ 10 15 cm 
and Vi are the gas-phase number density and thermal velocity of species i. 
Solving for the dust temperature we can derive the sublimation temperature 



for species i (Hollcnbach et al. 2008) 



Eb.i 
~k~ 



57 + In 



N s 



10 15 cm- 



10 13 s" 1 



1 cm 



(35) 

The above equation assumes that there is at least one monolayer of the par- 
ticular species present on the grain surface. In the case where only a partial 
monolayer exists then the thermal desorption rate per unit grain surface area 
is f s ,iN Si ik evap ,i, where / s ,j is the fraction of the N s ira 2 surface sites occupied 
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by species i. In Fig. ^b) we show a plot of the CO evaporation timescale 
as a function of grain temperature. Comparing to the freeze-out timescale in 
Fig. [7ja) suggests that CO will be frozen in the form of ice throughout much 
of the outer (r > 20 — 40 AU) disk midplane. However, on the warmer disk 



surface (see Fig. 12 1 CO will remain in the gas phase. This is consistent with 
observations ( |van Zadelhoff et al. 20011. 

Another interesting aspect of Eqn. |35|is the dependence on gas density. In 
Fig. [8] we provide a plot of the sublimation temperature (and dust tempera- 
ture) as a function of disk pressure and radius for a model solar analog disk 



(taken from D'Alessio et al. 2001 1. This calculation uses the binding energies 



given in Table |3j Thus most volatile ices are frozen beyond tens of AU with 
various species returning to the gas in a fashion that depends on their relative 
binding strength to the grain surface. It should also be noted that experi- 
ments suggest that the return of species to the gas is not entirely as simple 
as described above. The evaporation also depends on the composition of the 



mantle and the relative disposition of its components (Viti et al. 2004) 



Photodesorption : 

The direct absorption of an ultraviolet (UV) photon by a molecule ad- 
sorbed on the grain surface puts an electron in an excited state. If the inter- 
action between the excited state and the binding surface is repulsive, then the 



excited molecule may be ejected from the surface (Watson 1976 1. The rate of 



photodesorption per adsorbed molecule is given by (Boland & de Jong 1982 



Hollenbach et al. 2008 1 



I? 



pd 



Tra 2 d eYG F exp(-l.8A v ) cm 



-v 1 . 



(36) 

In this equation Y is the photodesorption yield (number of molecules desorbed 
per incident photon) , e is the fraction of species that are found in the top few 
monolayers, is the grain radius, and Go is the UV field enhancement factor 
in units of the mean interstellar UV field of Fq — 10 s photons cm -2 s _1 



(Habing 1968|). The stellar UV field and its relation the interstellar radiation 

For most molecules the yield is unknown and is 
3 (Bergin et al. 19951. More recently, there have 



4.2 



io- 



field is discussed in 
assumed to be ~ 10 
been measurements for CO (Y - 
and H 2 (V = T/100 x 1.3 x 10 



3 x 10~ a mol/photon; Oberg et al. 2007) 



mol/photon; Oberg et al. 20091 



X-ray Desorption: 

Young stars are known to be prolific emitters of X-ray photons with emis- 



sion orders of magnitude above typical levels for main sequence stars (Fcigcl 



son & Montmerle 1999). This emission is capable of penetrating much deeper 



into the disk atmosphere than UV photons and hence can potentially pro- 
vide energy to release molecular ices from grains on the disk surface. For this 
chapter we will focus on the effects of these X-rays on the chemistry and re- 
fer the reader to the chapter by Feigelson for a discussion of the origin and 
characteristics of X-ray emission. X-ray photons are absorbed by the atoms 
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Fig. 8. Plot of the sublimation temperature for various molecular ices, and the 
midplane dust temperature, as a function of disk pressure and radius for a disk 



model taken from D'Alessio et al. (2001) 



that compose the grain refractory core and mantle with a cross-section that 
depends on species and wavelength ( Morrison & McCammon 1983| |Henke| 



et al. 19931. The absorption of an X-ray produces a hot primary electron 



that generates a number of secondary electrons (s) via a process known as the 
Auger effect. Depending on the energy of these electrons, and the grain size, 



they then deposit all or a fraction of their total energy to the grain ( Dwek & 



Smith 1996 1. This heat diffuses through the grain lattice based on the thermal 



physics of the material. Fortunately experiments have demonstrated that the 
thermal conductivity and specific heats of amorphous solids tend have near 
universal values to within an order of magnitude over 3 decades in tempera- 
ture ( |Pohl| [1998] ). This is encouraging for interstellar studies as we can more 
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reliably adopt laboratory estimates for studies of grain heating (see appendix 
in 



Najita et al. 2001 ) 



Leger et al. ( 1985 1 explored grain heating by X-rays and cosmic rays adopt- 



ing the specific heat measurements of amorphous silica from Zcll er fc Pohl| 
( 1971 1 in the following form: 



pC v (T) = 1.4 x 10~ 4 T 2 J cm- 3 K~ x (10 < T < 50 K) (37) 
pC v {T) = 2.2 x lO^T 1 ' 3 J cm- 3 K~ l (50 < T < 150 K) (38) 

Using this expression the energy increase due to an energy impulse raises the 
temperature by (using the specific heat appropriate for T < 50 K): 



AE(eV) = 3.4 x 10 4 



0.1/im 



T. 



f 



30 K 



where is the final grain temperature and a is the grain size (Leger et al. 



19851. Najita et al. (20011 explored the question of X-ray desorption in pro- 



toplanetary disks using the results of |Dwek fc Smith (19961 and Leger et al 



( 1985 1 . As one example they found that a 1 keV photon will deposit ~ 700eV 



in a 0.1 /xm grain. The warm grain cools via evaporation and thermal radia- 
tion. This energy impulse is insufficient to raise a cold (Td us t ~ 10 K) grain 
to above the threshold for desorption of the volatile CO molecule. However, 
smaller grains (a < 0.03^m) have reduced volume for heat diffusion and equi- 
librate to temperatures that evaporate CO molecules. Thus X-ray desorption 
from the whole grain is not feasible, as grains with sizes > 0.1 /im are present 
in the ISM and in disks. A more efficient method for desorption is found 
when one considers the fact that grains are likely porous fractal aggregates 
composed of a number of small grains with a range of sizes. In this case the 
absorption of an X-ray by a portion of the grain is capable of local evapora- 



tion prior to the heat diffusing throughout the aggregate (Najita et al. 2001 1 



This is efficient for desorbing the most volatile ices (e.g. CO, N2) on the disk 
surface, but not water ice. There is some evidence that this process is active 



( Greaves , 2005 1 ; although the situation is complex because one additional ef- 
fect of X-ray ionization is the excitation of H 2 which generates a local UV 
field that could release molecules via photodesorption. 

Cosmic ray Desorption: 

Similar to X-rays, cosmic rays absorbed by dust grains can provide non- 
thermal energy to desorb molecules frozen on grain surfaces. Cosmic rays 
also have greater penetrating power, when compared to X-rays and UV, and 
hence can have a greater impact on the chemistry in well-shielded regions. 
The penetration of cosmic rays to the inner tens of AU of the disk is highly 
uncertain; however, it is likely that cosmic rays can affect the chemistry of the 
outer disk. This issue will be discussed later in § |4.2| and here we will focus on 
the effects of non-thermal desorption via cosmic ray impact. Cosmic rays have 
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energies between 0.02 - 4 GeV per nucleon, but the energy deposition goes 
as Z 2 such that the less abundant heavier cosmic rays have greater impact 



on the grain thermal cycling and ice evaporation (Leger et al. 19851. Bringa 



& Johnson ( 2004 1 noted that the impact of the cosmic ray on a grain is not 



exactly a question of thermal diffusion (as normally treated) , but rather one 
of a pressure pulse and a melt. In general, the effect of cosmic rays is to 
efficiently remove all ices (including water) from small (0.01 /im) grains due 
to large thermal impulses. Cosmic ray impacts can also release ices from larger 
grains via local spot heating. This is thought to be capable to remove the most 



volatile ices (CO, N 2 ), but not water ice (Leger et al. 1985 Herbst Cuppcn 



2006 ) . A number of groups have looked at this process and there is disparity 



in the calculations. A summary of the various approximations is provided by 



Roberts et al. (20071 



Explosive Dcsorption: 

Experiments of ices irradiated in the lab suggest that when an ice mantle 
formed at ~ 10 K is exposed to radiation and then warmed to ~ 25 K , the 
rapid movement (§ 3.2 1 and reactivity of radicals leads to an explosion of 



chemical energy and desorption (Schutte & Greenbergl 19911. There are two 



issues with this mechanism that are not insurmountable, but make it more 
difficult to readily place into models. (1) The generation of radicals on the 
surface needs the ice-coated grain to be exposed to the UV radiation field. A 
criticism of this is that the UV flux in the lab was significantly greater than 
typical exposure in the well shielded ISM where ice mantles form. However, 
the strength of the UV (Ly a) field at 10 AU in TW Hya, th e closest young 
star with a circumstellar disk, is ~ 10 14 photons cm -2 s _1 (Herczeg et al. 



2002 ) , comparable to that used in the laboratory. The UV radiation may also 



come from external high mass stars if a young solar-type star is born in a 
cluster. Thus it is likely that - on the disk surface - that radicals are being 
created in the frozen mantle. (2) Once the radicals have been generated then 
a cold (~ 10 K) grain needs to be warmed to above the temperature of ~25 K. 
This could be possible if the grain absorbs an X-ray, cosmic ray, undergoes a 
grain-grain collision, or advects to warmer layers. 



Grain Surface Chemistry 



A full treatment of grain surface chemistry is beyond the scope of this chapter 
and the reader is referred to the review bylHerbst et al 



(2005) 



There are a 

few issues specific to disk chemistry that should be discussed. First at low 
dust temperatures (< 20 K) H atom addition will dominate as H (and D) 
atoms can more rapidly scan the surface when compared to heavier atoms or 



molecules. The scanning timescale is (Tielens 2005): 



-E m /kT dv 



(39) 
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v m is the vibration frequency of the migrating molecules which is equivalent 
to Uq (Eqn.|33]). It is typical to assume E m ~ 0.3 — 0.5E b , thus at T dust ~ 10 K 
atoms with low binding energies (H, D) migrate rapidly and heaver molecules 
tend to remain in place. 

The H atom abundance in molecular gas is set by the balance between 
H2 ionization (which ultimately produces a few H atoms per ionization) and 
reformation of H2 on grains. Since the ionization rate in the disk midplane is 



likely small (see § 4.2) and cold grains are present, there may be few H atoms 



(or other atoms) present in the gas to allow for an active surface chemistry. 
To be more explicit, the sequence of H formation and destruction reactions in 
a well shielded molecular medium arc as follows: 



H + gr^H gr (40) 
Ccr,xr,rn+H 2 -> H+ + e (41) 
H+ + H 2 -> H+ + H (42) 

The limiting step in H atom formation is the ionization from cosmic rays 
(CR), X-rays (XR), or radionuclides (RN). In essence every ionization of H 2 
ultimately produces fn ~ 2-3 H atoms per ionization. In equilibrium the H 
atom concentration is: 

n H= 5 — = Tn=2i c — ■ ( 43 ) 

n g OgV H aHr\ 10 ll V H S H T] 

In this equation Vh is the thermal velocity of hydrogen atoms, and 77 is an 
efficiency factor which can be set to unity provided that there is more than 
one H atom on each grain. Sh is the sticking coefficient for H atoms. |Burke fc] 
HoUenbach ( |1983[ ) and |Buch fc Zhang| ( |1991| ) provide analytical expressions 



that can be used with values of Sh ~ 0.7 at 20 K. 



We have also used the expression for n g a g = 10 n discussed in § 3.2 



Each of the relations in Eqn. 43 are constant (given a temperature), thus the 
space density of H atoms is a constant in the molecular gas. If we assume that 
cosmic rays penetrate to the midplane with an ionization rate of £ = 1 x 10 - 
s _1 (the most extreme high ionization case for the midplane) and the gas 
temperature is 30 K then nn ~ 5 cm~ 3 . The H atom abundance is therefore 
inversely proportional to the density. The average grain in the ISM with a 
size of 0.1 fxm has an abundance of ~ 10~ 12 . This sets the initial conditions 
for the abundance of grains. If 77 = 1 then for much of the disk midplane, 
where n > 10 12 cm -3 there will be less than 1 H atom per grain significantly 
reducing H atom surface chemistry. However, it may be the case that there is 
less than one H atom on each grain and 77 <gc 1 which could significant reduce 
the grain surface formation rate. This would lead to more gas phase H atoms 
and a potential to power chemistry but requires longer timescales. 

Second, at higher temperatures heavier atoms and molecular radicals can 
begin to diffuse on the grain surface and create complex species. Given the 
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abundance of radiation on the disk surface it is highly likely that bonds are 
broken in the ices creating radicals (from CH 3 OH, H 2 0, CO2, CH 4 , NH 3 , ...). 



Upon warm up this can produce more complex species (see models of Garrod 



et al. 2008 as an example), unless it is limited by explosive desorption. Thus 
grain chemistry may be important on the disk surface and, depending on the 
strength of vertical mixing, its products may reach the disk midplane. 



Deuterium Fractionation 

Enhancements of deuterium-bearing molecules relative to the hydrogen coun- 
terparts have been known for some time in solar system and in the interstellar 
medium. In cold (Tdust ~ 10-20 K) cores of clouds where stars are born, en- 
richments of 2-3 orders of magnitude are observed above the atomic hydrogen 
value of (D/H) > (2.3±0.2) x 10" 5 estimated within 1 kpc of the Sun (Linsky 



et al. 2006 ) . Because of the lower zero point energy for deuterium compared 



to hydrogen, at low (T < 30 K) temperatures, it becomes favorable to transfer 
the D-bond, as opposed to the H bond. Ion-molecule reactions in the dense 



ISM are thought to be the mechanism responsible for these enrichments (Mil 



lar et al. 1989). Deuterium fractionation can also be powered via reactions 



2005), which is 



on the surfaces of dust grains (Tielens 1983 Nagaoka et al. 
also discussed below. 

In the gas, deuterium chemistry is driven by the following reaction: 



Ho 



HD <-> H 2 D+ + H 2 + 230K. 



(44) 



The forward reaction is slightly exothermic favoring the production of H 2 D + 
at 10 K, enriching t he [D]/[H] rat io in the species that lie at the heart of ion- 
molecule chemistry ( Millar et al.[ [l989p| These enrichments are then passed 
forward along reaction chains to DCO + , DCN, HDO,HDCO, and others. 

At densities typical of the dense interstellar medium (n ~ 10 5 cm -3 ), pure 
gas-phase models without freeze-out cannot produce significant quantities of 
doubly (NH D 2 ; |Roueff et al. 2000[) a nd triply deuterated ammonia (ND 3 ; Lis 
et al.| 2002| van der Tak et al. "[2002). This motivated a re-examination of the 
basic deuterium chemistry. The primary advance in our understanding is two- 
fold: (1) deuteration reactions do not stop with H 2 D + , rather they continue 
towards the formation of both D 2 H + and D3", via a similar reaction sequence 
(Phillips & Vastell [2003] IRoberts et al.l [20031 IWalmsley et all [20041): 



H 2 D+ + HD ^ D 2 H+ + H 2 + 180K, 



(45) 



H 2 D+ + HD ^ D+ + H 2 + 230K. 



(46) 



4 Reaction 44 has been measured in the lab at low temperatures finding that there 
is an additional dependence on the ortho/para ratio of H2 (Gerlich et al. 20021. 
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(2) The freeze-out of heavy species in the disk midplane (§ 3.2 1, in particular 



CO, a primary destroyer of both Hjj~ and H 2 D + , increases the rate of the gas 
phase fractionation reactions. Thus in vertical gas layers, prior to potential 
complete freeze-out of heavy molecules, there can exist active deuterium frac- 



tionation, provided the gas temperature is below 30 K (see Aikawa & Herbst 



2001) 



At high densities in the inner nebula thermodynamic equilibrium can ap- 
ply. At temperatures T > 500 K [HDO]/[H 2 0] closely follows [HD]/[ H 2 ], but 
at T < 500 K fractionation can proceed via the following reaction (Richet 



et al. 1977): 



HD + H 2 ^ HDO + H 2 



(47) 



This reaction can provide a maximum enrichment of a factor of 3 (Lccluse & 



Robert 19921 



Grain surface chemistry can also produce deuterium fractionation as the 
gas-phase chemistry (at low temperature) leads to an imbalance and enhance- 
ment in the D/H ratio of atomic H. The efficiency of this process depends on 
the number of atoms which arc generated by ionization. At low temperatures 
more deuterium is placed in H 2 D + relative to H^~. These species create H and 



D atoms by reacting with neutrals. Thus fr> > fn (see Eqn. 43 1 and there 



is the potential for these atoms to freeze onto grains and react, thereby frac- 
tionating molecular ices. However, as noted in § |3.2| it is possible that in the 
dense midplane that there will be less than one D atom generated per grain 
and the efficiency of D-fractionation on grain surfaces could be reduced. 

To summarize, deuterium chemistry can be active in the disk, both in the 
gas and on grains, but it requires low temperatures and sufficient levels of 
ionization. 



3.3 Kinetics and Thermodynamic Equilibrium 



At present it is not clear exactly where in the disk thermodynamic equilibrium 
is valid. It is clear that some aspects of meteoritic composition are consistent 
with equilibrium condensation (§ 2.1 1. However, there are also aspects that 



are inconsistent, such as the detection of unaltered pro-solar grains (Huss 
|fc Lewis) 1995 1 and deuterium enrichments in cometary ices, which suggest 
that not all material cycled to high temperatures (above evaporation). Fegley 
( 2000 1 provides a schematic of the relative importance of equilibrium to Id 



netics. He proposes that thermodynamic equilibrium is of greater significance 
in the inner nebula and within giant planet subncbula. To some extent this 
must be the case, but divining the exact transition is not entirely clear as, 
for example, the inner nebula and surfaces of giant planet subnebula will be 
exposed to radiation that powers non-equilibrium photochemistry. Given ev- 
idence from trends in meteoritic composition (e.g. Davis & Richter 2004 as 



one example) , it is clear that in some instances thermochemical equilibrium is 
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valid. However, one must pay careful attention, where possible, to the relevant 
kinetic timescales. 



4 Physical Picture and Key Processes 



4.1 Density and Temperature Structure 

Estimates of the physical structure of circumstellar disks can come from the 
solar system and also from extra-solar systems. Within the solar nebula esti- 
mates are based on the so-called minimum mass solar nebula: the minimum 
amount of mass required to provide the current distribution of planetary mass 



as a function of distance from the Sun, correcting for a gas/dust ratio. Hayashi 
( 1981 1 estimates E(r) ~ 1700r -3 - 2 , which is commonly used, but see also Wei- 



denschilling (19771 and Davis (20051. It is important to note that this value 



is a lower limit and does not account for any potential solids that have been 
lost from the system during the early evolutionary phases. Limits on the gas 



temperature have been estimated by Lewis (19741, based on the differences in 
planetary composition. For instance the temperature at 1 AU is suggested to 
be ~ 500 — 700 K, based on the estimation that the Earth formed below the 
condensation temperature of FeS, but above the end point for oxidation of 

Based on these estimates the nebular 
which is much steeper than in pure 



metallic iron in silicates (|Lewis| 1974). Based on these estimates the nebular 

-l 



thermal structure would decay as R 
radiative equilibrium R~ 5 (see, e.g. Hayashi 1981). However, as discussed 



below, radiation does dominate the thermal structure of protoplanetary disks. 

More detailed estimates as a function of position (both radial and verti- 
cal) are available from extra-solar protoplanetary systems. Models take into 
account the energy generated by accretion and by stellar irradiation of a 



flared disk structure. Disk flaring is due to hydrostatic equilibrium (Kcnyon 



& Hartmann 1987) and is characterized by a pressure scale height H p 



^/R 3 c^/GM^, with the sound speed, c s = y/kT/pmn (p is the mean molec- 
ular weight, p ~ 2.3). The vertical density then follows p = poexp(— z 2 /2H 2 ) , 
with po a function of radius. With radial density profile as a variable and 
some assumptions regarding dust composition and size distribution (to set 
the optical properties) the spectral energy distribution of dust emission from 
UV/Visible (~ 4000 A) to millimeter (1000 pm) can be predicted and com 
pared to observations ( Calvet et al. 19911 Chiang & Goldreich 1997 D'Alessio 



et al. 



19981. These models have found that disk surfaces, which are directly 
are warmer than the midplane. The midplanc 
for the inner few AU, by vis- 



heated by stellar irradiation 
is heated by reprocessed stellar radiation and 
cous dissipation of accretion energy. Fig. [9] shows an example of the thermal 



structure from one such model (D'Alessio et al. 2001) 



An important issue is the effect the initial steps of planetesimal forma- 
tion, grain growth and settling, on the disk physical structure. Smaller grains 
((a) ~ 0.1 pm) preferentially absorb the energetic short wavelength and as 
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Fig. 9. Predicted te mperature structure for a typical T Tauri disk model by 



D'Alessio et al 



(2001 1. Assumed parameters are 0.5 Mq, M = 3 x 10" s M Q /yr, and 



a = 0.01. 



these grains grow to larger sizes their optical depth decreases (Dullcmond & 



Dominik 2004 D'Alessio et al. 20061. The presence of small grains in the 



disk atmosphere leads to higher temperatures and greater flaring than disks 
with larger grains (or more settling to the midplane). Thus the disk physical 
structure evolves with the evolution of solids. Moreover there is an observed 
decay in accretion rate over similar timescales (iMuzerolle et al. 20001. Thus, 



the energy input from accretion should similarly decline. In sum, both the 



midplane temperature and pressure should decline with evolution (D'Alessio 
et al. 2005 ) . As we will show in § [5]this structure and evolution is an important 
component of time-dependent disk chemical models. 
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4.2 Disk Ionization 



The gas-phase chemistry is powered by ionizing photons. The characteristics 
of the chemical kinetics therefore depend on the flux and wavelength depen- 
dence of ionizing photons. The molecular ions generated by photo- or cosmic 
ray absorption are also the key to linking the mostly (by many orders of mag- 
nitude) neutral gas to the magnetic field which is thought to have a direct 
relation to the physics of disk accretion. In the following we will explore each 
of the various components of disk ionization. In Table [4~2| we provide a listing 



of the major contributors to the ion fraction (table taken from Bergin et al. 
20071. 



Cosmic Rays 

The interstellar cosmic ray ionization rate has been directly constrained by 
the Voyager and Pioneer spacecraft at ~ 60AU from the Sun ( Webber , 1998 1 



and estimates in the ISM have been summarized by Dalgarno ( 2006 1 . Based 
on these studies, the ionization rate is believe to be (jy 2 ~5x 10 1Y s _1 in the 
dense ISM with growing evidence for higher rates in the diffuse gas (McCall 



et al. 2003 Indriolo et al. 2007 1 . It is this rate that impinges on the surface of 



the disk. Umebayashi & Nakano ( 1981 ) explored penetration depth of ionizing 



protons as a function of gas colu mn a nd find that cosmic rays penetrate to 
a depth of S ~ 96 g cm -2 . Table 4.2 compares this e-folding depth to that 



of other potential contributors to the ion fraction. Excluding radionuclides, 
cosmic rays represent the best mechanism to ionize and power chemistry in 
the disk midplanc. However, if present, even cosmic rays will have difficulty 
penetrating to the midplanc in the inner several All's of disks around solar- 
type stars. Another important question is whether ionizing cosmic rays are 
present in the inner tens of AU at all. Within our own planetary system, the 
Solar wind limits ionizing cosmic rays to beyond the planet formation zone. 
Estimates of mass loss rates from young star winds significantly exceed the 



Solar mass loss rate (Dupree et al. 20051, and may similarly exclude high 



energy nuclei. Cosmic rays will contribute to the chemistry for the outer disk 



and the detection of ions potentially provides some confirmation (Ceccarelli 



& Dominik 2005 ) 



Ultraviolet Radiation 



T Tauri stars are known to have excess UV flux much higher than their ef- 
fective temperature, T e // ~ 3000 K (Herbig & Goodrich 19861, which are 
generated, at least in part, by accretion ( Calvet fc Gullbring||1998 1. Observa- 
tions have suggested that the UV radiation field (and X-ray ionization) play a 
key role in the observed molecular emission ( Willacy fc Langerj |2000| Aikawa 



et al. 2002 1. UV radiation is important as the primary molecular photodisso- 



ciation bands lie below 2000 A and this radiation both dissociates molecules 
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and can power a rich chemistry on the disk surface (sec, for example Ticlcns 



& Hollenbach 19851 



In this context there are 2 photon fields to consider. The stellar radiation 
field and the external radiation field which can be enhanced when compared 
to the standard interstellar radiation field (ISRF). It is common to place both 
the stellar and interstellar UV radiation field in the context of the ISRF which 
is defined by the measurements of Habing (1968) and Draine (1978). Using 
the measurements of Habing, the energy density of the ISRF is 1.6 x 1CP 3 
erg cm -2 s _1 which is defined as Go = 1 (in this context the Draine field is 
G = 1.7). As an example the FUV flux below 2000 A for T Tau is - 3 x 10" 13 
erg cm -2 s" 1 A -1 . Scaled to 100 AU at 140 pc and integrated over 1000 A (the 
Lyman limit to 2000 A where most molecular photoabsorption cross-sections 
lie) gives 23.5 erg cm" 2 s" 1 or G = 1.5 x 10 4 at 100 AU (scaling as 1/r 2 ). 
T Tau has a particularly high mass accretion rate of dM/dt ~ 10~ 7 M /yr 
that is about an order of magnitude higher than typically observed. Most T 
Tauri stars have lower accretion rates and weaker UV fields with values in the 
range of G o (100 AU) = 3 - 1000, based on observations from FUSE and 
HST ( Bergin et al.| 2004). Bergin et al. (20031 pointed out the importance 
of Ly a emission, which is present in the stellar radiation field, but absent in 
the ISRF. In one star that is unobscured by interstellar absorption, TW Hya, 



Ly a contains as much of as 85% of the stellar UV flux (Herczeg et al. 



2002 1. Even in systems where Ly a is absent due to absorption and scattering 



by interstellar H I, Ly a pumped lines of H 2 are detected, testifying to the 



presence of Ly a photons in the molecular gas (Herczeg et al. 2002 Bergin 



et al. 


2004 


Herczeg 


et al. 


2006) 


stars is shown in Fij 


10 





A sample of the UV field of typical T Tauri 



Stellar photons irradiate the flared disk with a shallow angle of incidence, 
while external photons impinge on the disk at a more normal angle. Inter- 



stellar photons therefore have greater penetrating power (Willacy & Langer 



20001. van Zadelhoff et al. (2003) employed a 2D model of UV continuum 



photon transfer and demonstrated that scattering of stellar photons is key to 
understanding how UV radiation influences disk chemistry. Fig. 11 illustrates 
the geometry of various contributors to the radiation field on the disk surface. 
An analytical approximation of the transfer of the stellar field is provided by 
Bergin et al. ( 2003 1 . In Taurus the average extinction towards T Tauri stars 



is Ay ~ l m ( Bergin et al.[ 2003); under these circumstances the stellar field 
will dominate over the ISRF (Go = 1) for much of the disk. 

Using the parameterization of the stellar field with respect to the ISRF 
allows the use of photodissociation rates calculated using the strength and 
shape of the IRSF, assuming albedo and scattering properties estimated for 
interstellar grains (van Dishoeck 1988 Roberge et al.[[l991 l. This correctly 
treats the direct attenuation of the field, but not the contribution from scatter- 
ing, which dominates the deep surface (van Zadelhoff et al. 2003 1. The shape 
of the stellar field, which is primarily line emission, is also different from the 
ISRF (see Fig. 10 1. This is particularly acute for Ly a as some species are 
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Fig. 10. UV spectra of T Tauri stars. Heavy solid lines and light solid lines represent 
the spectra of BP Tau and TW Hya, respectively. The spectrum of TW Hya is scaled 
by 3.5 to match the BP Tau continuum level. The long dashed line represents the 
interstellar radiation field of Draine (1978) scaled by a factor of 540. The region 
around the Lya line is enlarged in the right panel. Taken from Bergin et al. (2003). 



sensitive to dissociation by Ly a emission (e.g. H 2 0, HCN), while others are 
not, such as CO and CN ( Bergin et al.[|2003| van Dishoeck et al. 



2006). In 



this case, it is better to directly calculate the photodissociation and ionization 
rates from the flux and the cross-sections. This can be done using the observed 
UV field, or a proxy black-body, in the following fashion: 



^photo 



4ttA 
he 



a(X)J x>0 (r)e-^dX (s" 1 ). 



(48) 



er(A) is the photodissociation cross-section and Ja,o is the mean intensity of 
the radiation field at the surface (Fx = 4-7T J\ = J I\df2, where i? is the solid 
angle of the source). The main source of opacity at these wavelengths is dust 
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External radiation: - 
normal angle of incidence 



Fig. 11. Schematic showing the difference in angle of incidence between stellar 
and external radiation impinging on a flared disk. Radiation with a shallow angle 
of incidence will have less penetration into the disk. Scattering of stellar radiation 
increases its importance, particularly if the external radiation field suffers from any 
local extinction from a surrounding molecular cloud. 



grains and some assumptions must be made regarding the composition, opti- 
cal properties, and size distribution (see §4.3 for a more complete discussion 
of dust grain absorption of radiation), van Dishoeck et al. (20061 provides 
photo-rates for various blackbodies and also photodissociation cross-sections 
coincident with Ly a. In addition, a web site for cross-sections is available at 
http : / / www . strw . leidenuniv . nl/~ewine/photo] 

It is important to note that the photodissociation of two key species, H2 
and CO, require additional treatment. These species are dissociated via a line- 
mediated process, as opposed to a dissociation continuum typical for most 
other molecules. In this regard the optical depth of the lines depends on 
the total column. When the lines become opaque molecules closer to the ra- 
diation source can shield molecules downstream in a process that is called 
self-shielding. In addition, H2 molecules can partially shield some of the CO 
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predissociation lines. Self-shielding rates are time and depth dependent. In 
general, approximations based on the interstellar UV field are adopted from 



Draine & Bertoldi (1996), for H2, while for CO the shielding functions pro- 



vided by van Dishoeck & Black (19881 and Lee et al. (19961 can be used. 
One caveat of these approximations is that the shielding functions are calcu- 
lated assuming that molecules predominantly reside in the ground state with 
a line profile appropriate for cold interstellar conditions. In regions with high 
temperatures these approximations can potentially break down. 



X-ray Ionization 



X-rays are a key source of ionization and heating of protoplanetary disks. The 
typical X-ray luminosity (based on protostars in Orion) is L x = 10 28 5 — 10 31 
ergs s _1 with a characteristic temperature of the X-ray emitting plasma of 
T x ~ 1 — 2 keV (e.g. Feigelson & Montmerle 1999 see also chapter by Feigel- 



son, this volume). X-ray photons interact directly with the atoms in the disk; 
i.e., the inclusion of the atom in a molecule does not affect the photoabsorption 
cross-section. In particular, the opacity is dominated by inner shell ionization 
of heavier elements. For atoms heaver than Li the inner shell ionization is 
followed by the Auger effect which generates a number of primary electrons. 
Each primary electron (p) produces secondary electrons (s) by impact ion- 



ization of the gas, with N s 



^> N p . Therefore for the X-ray ionization rate, 



C* = C P + Cs = C- (49) 
This is a well known result and is similar to that seen for cosmic rays. In 



general there are N s ~ 30 atoms and molecules ionized per keV (Aikawa & 



Hcrbst 2001 1 



The X-ray ionization rate can be derived by an integration of the X-ray 
flux, F(E), and the cross-section tJi(E) of element i (in this case hydrogen) 
over energy: 



30 keV 



- N s / ai(E)F{E)dE (50) 

J l keV 

The X-ray flux can be estimated from the X-ray luminosity assuming a charac- 



teristic temperature for the plasma (see Glassgold et al. 1997 for an example) . 
The X-ray flux decays with depth into the disk, F(E) = Fxfie~ T ^ E \ where 
t(E) = Nctt{E). X-ray photons are attenuated by photoabsorption at the 
atomic scale. In this regard is customary to use the total X-ray cross-section, 
<7t(E), computed using solar abundances (Morrison & McCammon 19831. 



With this assumption the total hydrogen column is the only information re- 
quired. Table 5] list the fit to the X-ray cross-section computed by |Morrison fc| 
McCammon (19831. This assumes that the heavy elements, that control the 



opacity, have Solar abundances and are uniformly distributed in the disk. 
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Table 5. Coefficients for Fit to X-ray 
Cross-Section a ' b 



Energ 


y Range 


(kcV) c 




Cl 




C2 


0.030 


-0.100 


17 


3 


608 


1 


-2150.0 


0.100 


-0.284 


34 


6 


267 


9 


-476.1 


0.284 


-0.400 


78 


1 


18 


8 


4.3 


0.400 


-0.532 


71 


4 


66 


8 


-51.4 


0.532 


-0.707 


95 


5 


145 


8 


-61.1 


0.707 


-0.867 


308 


9 


-380 


6 


294.0 


0.867 


-1.303 


120 


6 


169 


3 


-47.7 


1.303 


-1.840 


141 


3 


146 


8 


-31.5 


1.840 


-2.471 


202 


7 


104 


7 


-17.0 


2.471 


-3.210 


342 


7 


18 


7 


0.0 


3.210 


-4.038 


352 


2 


18 


7 


0.0 


4.038 


-7.111 


433 


9 


-2 


4 


0.75 


7.111 


-8.331 


629 





30 


9 


0.0 


8.331 


-10.000 


701 


2 


25 


2 


0.0 



Fit and tabulation reprinted from |Mor^] 



rison fc McCammon (1983) 
b Cross-section per hydrogen atom, 
<tt{E), is (co + dE + aE'^E- 3 x 10~ 24 
cm 2 (E in keV). 



It is well known that dust grains suspended in the disk atmosphere will 
sink to the midplane, reducing the dust/gas mass ratio in the upper disk at- 
mosphere ( Weidenschilling & Cuzzi 1993). This will reduce the optical depth 
of the upper layers to both X-rays and UV radiation. Compared to UV trans- 
fer, dust coagulation will not reduce the X-ray opacity as dramatically, unless 
grains grow to large sizes. Instead the mass of heavy elements must be redis- 
tributed as is the case of dust settling. In the extreme limit only H and He are 
present in the upper atmosphere and the cross-section at 1 keV is reduced by 
a factor of ~4.5 (i.e. the H and He absorption limit, Morrison & McCammon 



19831. In the case of non-solar abundances or mass redistribution the cross- 
section for each individual element can be used (Henke et al. 19931. X-ray 



ionization is crucial for the inner disk where cosmic rays likely do not pene- 
trate. X-ray photons also have a significantly greater penetration depth than 



UV photons and thus can power chemistry (and perhaps accretion, Gammie 



1996 1 on the disk surface. The characteristics of the X-ray driven chemistry 



can extend beyond the simple ionization of H2 as X-rays can produce short- 
lived, but reactive, doubly ionized molecules. For greater detail, the reader 



is referred to Langer (19781, Krolik & Kallman (1983), and Stauber et al. 
(20051. 
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Active Radionuclides 



In the case that the midplane of the disk is completely shielded from cosmic 
ray radiation then radionuclides provide a baseline level of ionization. A pri- 
mary issue in this regard is the overall abundance of various radionuclides 
and whether these elements are uniformly distributed throughout the sys- 
tem. Key elements in this regard are 26 A1 and 60 Fe which have short lifetimes 
[t 1/2 ( 26 A1) = 0.72 Myr; T 1/2 ( 60 Fe) = 1.5 Myr ] and have been inferred to b e 
present in the solar system (e.g. Wadhwa et ah] 2007 Wasserburg et al. 20061. 
If present they are a strong source of ionization for several millions of years. 
Table |4.2| provides an estimate of this ionization rate, which is dominated by 
26 Al. This is sufficient to allow for an active ion driven gas chemistry in regions 
where neutral molecules (excluding H2 and HD) are not completely frozen on 
grains. This ionization rate will decay with time according to respective half- 
life. If 60 Fe or 26 Al are not present then this rate gets substantially reduced. 



Finocchi & Gail (1997) provide a short summary of the relevant ionization 



processes including the most important radionuclides. 



Dissociative Recombination 



A key to the question of the overall ionization fraction is the balance between 
the ionization rate and the recombination rate of atoms and molecules. In equi- 
librium, the ion fraction, x e = n e /n, can be expressed by: x e = \J ( / (a r n) , 
where a r is the electron recombination rate. The dependence of the ion frac- 
tion on position, therefore not only depends on the flux of ionizing agents, but 
also on the recombination rate of the most abundant ions. There are some key 



differences in this regard that are summarized in Table 4.2 For the midplane 
a key question is whether molecular ions, metal ions, or grains are the domi- 
nant charge carriers and this serves as a useful foil to discuss some differences 
in recombination. 

In general atomic species have longer recombination timescales than molec- 
ular ions. Thus the presence of abundant slowly recombining metals can have 
key consequences for magnetic field coupling ( Ilgner & Nelson 2006 1 . How- 



ever, metal ions are not present in the dense star forming cores that set the 



initial conditions for disk formation ( Maret & Bergin 2007 ) , and are likely not 



present in the disk gas. If grains are the dominant charge carrier the recombi- 
nation rate, a gr , is the grain collisional timescale with a correction for long- 
distance Coulomb focusing: a gr = ira 2 d n gr v(\ + e 2 / kadTdust) ■ At Tdust = 20 K, 



Draine & Sutin (1987) show that for molecular ions, grain recombination will 
dominate when n e /n < 10~ 7 (a m i n /3A) -3 / 2 . Grains can be positive or neg- 
ative and carry multiple charge: |Sano et"aL (20001 find that the total grain 
charge is typically negative, while the amount of charge is l-2e~, varying with 



radial and vertical distance. Florescu-Mitchell & Mitchell (2006) provides a 



summary of molecular ion recombination rates determined in the laboratory. 
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4.3 Grain Growth and Settling 



The onset of grain evolution within a protoplanetary disk consists of collisional 
growth of sub-micron sized particles into larger grains; the process continues 
until the larger grains decouple from the gas and settle to an increasingly 
dust-rich midplane (Nakagawa et al. 1981 Weidenschilling & Cuzzi 1993| ). 



Grain coagulation can alter the chemistry through the reduction in the 
total geometrical cross-section, lowering the adsorption rate and the Coulomb 
force for ion-electron grain recombination. Micron-sized grains couple to the 



smallest scales of turbulence (Weidenschilling & Cuzzi 19931 and have a ther- 



mal, Brownian, velocity distribution. Thus, the timescale of grain-grain colli- 



sions is, t. 



5/2 
- gr oc a j 



£ the gas-to-dust mass ratio, and T^ ust the 



dust temperature ( |Aikawa et al. 1999). In this fashion grain coagulation pro- 
ceeds faster at small radii where the temperatures and densities are higher. 
Aikawa et al.| ( |1999[ ) note that the longer timescale for adsorption on larger 



grains leaves more time for gas-phase reactions to drive toward a steady-state 
solution; this involves more carbon trapped in CO as opposed to other more 
complex species. 

Overall, the evolution of grains, both coagulation and sedimentation, can 
be a controlling factor for the chemistry. As grains grow and settle to the 
midplane, the UV opacity, which is dominated by small grains, decreases, 
allowing greater penetration of ionizing/ dissociating photons. As an example, 
in the coagulation models of iDullemond & Dominik ( 2004 ) the integrated 



vertical UV optical depth at 1 AU decreases over several orders of magnitude, 



towards optically thin over the entire column (see also Weidenschilling 19971. 

This can be understood by exploring the question of dust extinction of 
starlight from the perspective of the interstellar medium. In this case the 
amount of dust absorption is treated in terms of magnitudes of extinction, 
labelled as A^^At any given point in the disk (or in the ISM) the intensity 
can be charactered by the intensity (ergs s _1 cm~ 2 A -1 sr~ x ) that impinges 
on the surface, I V fi, modified by the dust absorption, I v = /„.oexp(— t>). The 
opacity, t\, is given by t\ — Nd ust Q\<j, where N^ust is total dust column 
along the line of sight, Q\ is the extinction efficiency, and a is the geometrical 
cross-section of a single grain. The wavelength dependent extinction, placed 
in terms of magnitudes, is defined below. 



A x = -2.b\og{I v /I v , Q ) (51) 

= 2.51og(e)r A (52) 

= \m&N dust Qxa (53) 

= 1.086 — NQ x a (54) 



Magnitudes in astronomy are defined such that a difference of 5 mag between 
2 objects corresponds to a a factor of 100 in the ratio of the fluxes. Thus an 
extinction of 1 mag corresponds to a flux decrease of (100) 1//5 ~ 2.512. 
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The total dust column can be related to the gas column (N) via the gas to 
dust mass ratio, f, which is measured to be 1% for interstellar grains (i.e. the 
primordial condition): 



Pdust 
Pgas 

n 



n d m gr 



nfirriH 

3 £,p,m H 
4ir p gr a 3 gr ' 



(55) 



(56) 



Placing the dust abundance (n c i us t/ n ) m t° Eqn. 54 and exploring the extinc- 
tion at visible wavelengths where Qy ~ 1, assuming a typical grain size of 0.1 
/itm, and p gr = 2 g cm -3 then, 



A v = 1.086 



3 £/j,m H 



0.01 



Qv 
1 



10" 21 iV. 



(57) 



Thus a total gas column of 10 cm provides 1 mag of extinction at visible 
wavelengths, for a standard gas-to-dust mass ratio. As grains settle to the 
midplane the gas remains suspended in the upper layers and the gas-to-dust 
ratio decreases. This requires larger columns for the dust to absorb stellar 
radiation. In our example we have explored extinction in the visual, while 
most molecules are dissociated by ultraviolet radiation. Dust extinction is 
larger in the UV (requiring a smaller gas column); however, the dependency 
on the gas-to-dust ratio and grain size remains. 

As grains evolve there will be a gradual shifting of the warm molecular 



layer deeper into the disk, eventually into the midplane ( Jonkheid et al. 2004) 



Because the grain emissivity, density, and temperature will also change, the 



chemical and emission characteristics of this layer may be altered ( Aikawa & 



Nomura 20061. These effects are magnified in the inner disk, where there is 



evidence for significant grain evolution in a few systems ( Furlan et al. 2005 1 



and deeper penetration of energetic radiation (Bergin et al. 20041. A key 



question in this regard is the number of small grains (e.g., PAHs) present 
in the atmosphere of the disk during times when significant coagulation and 
settling has occurred. 



4.4 Turbulence and Mixing 



It is becoming increasingly clear that mixing played an important role in the 
chemistry of at least the solids in the nebula. This is due to the detection 



of crystalline silicates in comets (e.g. Brownlee et al. 2006) and the presence 
of chondritic refractory inclusions in meteorites ( MacPherson et al.| |1988 1. 



Since the smallest solids (micron size grains) are tied to the gas it is likely 
that the gas is also affected to some extent (perhaps in a dominant fashion) by 
turbulent mixing. The question of mixing and its relevance has a long history 
in the discussions of solar nebula chemistry (e.g. Prinn 1993 and references 
therein) . 
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In terms of the dynamical movement of gas within a protoplanetary disk 
and its chemical effects, a key question is whether the chemical timescale, 
Tchem, is less than the relevant dynamical timescale, Td yn , in which case the 
chemistry will be in equilibrium and unaffected by the motion. If Td yn < Tchem 
then mixing will alter the anticipated composition. These two constraints are 



the equilibrium and disequilibrium regions (respectively) outlined in Prinn 
( |1993[ |. What is somewhat different in our current perspective is the recogni- 
tion of an active gas-phase chemistry on a photon-dominated surface (§ |5J. 
This provides another potential mixing reservoir in the vertical direction, as 
opposed to radial, which was the previous focus. 

It is common to parameterize the transfer of angular momentum in terms 
of the turbulent viscosity, v t — ac s H p , where v t is the viscosity, c s the sound 



speed, H p the disk scale height, and a is a dimensionless parameter ( Shakura & 



Sunyaev 1973 Pringle 1981 ). Hartmann et al. ( 1998 1 empirically constrained 
the a-parameter to be <; 10~ 2 for a sample of T Tauri disks. 

A number of authors have begun to explore the question of including 
dynamics into the chemistry. A brief and incomplete list is provided by |Bergin| 



et al. (20071. While the details differ, a common approach is to use mixing 



length theory where the transport is treated as a diffusive process (Allen et al. 
1981| jXie et al.[|1995| ), but see also |Ilgner et al.] ( |2004[ ); |Tscharnuter fc Gail 



(20071. For mixing length models the fluctuations of abundance of species 
can be determined by the product of the abundance fluctuations in a given 
direction and the mixing length (I), dxi ~ —Idxi/dz, where z denotes the 
direction where a gradient in abundance exists. The net transport is then 



(Willacy et al. 20061 



0i(cm 2 s 1 ) = nH 2 {v t dxi) 



Dun, 



dxj 
dz 



-Dm 



1 drii 
rij dz 



1 drift 2 



m 2 



dz 



(58) 

where D is the diffusion coefficient and v t is the turbulent velocity. The dif- 



fusivity is related to the viscosity v t by D = (vtl) 



ac s H p . Using this 



description, the chemical continuity equation can be written as: 



drij 
~dt 



d<f>j 
dz 



— Pi ~ Li 



(59) 



where Pi and Li are the chemical production and loss terms for species i. 



The radial disk viscous timescale is t v 



ioV 



a 



-l / t 

l gas 

UOOK 



' jv or, 
r 

1AU 



Mo 



The diffusivity, D or K, is not necessarily the same as the viscosity, v t (e.g. 
Stevenson 19901, as given above. Moreover, it is not entirely clear that ra- 



dial mixing will be the same as vertical mixing. In the case of disks, several 
groups have explored this question with a range of potential solutions ranging 
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from v t jD below unity to near 20 1 


Carballido et al. 


2005 


Johansen et al. 


2006 


Turner et al.| 2006 


Pavlyuchenkov & Dullemond 


2007 


|. For complete- 



ness, when Vt/D > 1 turbulent mixing is much less efficient than angular 
momentum transport. 



5 Current Understanding 

In this section we will synthesize the various physical and chemical processes 
in a discussion of our evolving understanding of disk chemistry. One important 



issue is that because of its low binding energy to grain surfaces (Hollcnbach 



& Salpeter 1971 1, molecular hydrogen will remain in the gas throughout the 
nebula. Provided that sufficient levels of ionizing agents reach the midplane or 
the deep interior of the disk surface then the ion-molecule chemistry outlined 
earlier can proceed. Surface chemistry based on diffusing hydrogen atoms 
will be more difficult to initiate, but in warmer regions heavier radicals can 
potentially migrate and react. Thus both gas and grain surface chemistry is 
important. 

For other molecules we expect large compositional gradients in the verti- 
cal and radial directions. The dominant effect is the sublimation/freeze-out 
as molecules transition from warm to colder regions in a medium with high 
densities [n ~ 10 7 — 10 15 cm~ 3 ) and therefore short collisional timescales be- 
tween gas and solid grains. In large part the overall chemical structure follows 
the thermal structure and is schematically shown in Fig. [12] This plot illus- 
trates that, beyond ~40 AU, the disk can be divided into three vertical layers. 
CO, as one of the most volatile and abundant species, controls the gas-phase 
chemistry in the outer disk and sets the radial and vertical boundary of these 
layers. The top of the disk is dominated by stellar UV and X-ray radiation 
which leads to molecular photodissociation in a photon-dominated layer. This 
layer will have several transitions. H/H2 is the first transition, as molecular 
hydrogen is the strongest at self-shielding. This is followed by C 11/C i/CO 
and subsequently the oxygen/nitrogen pools which require dust shielding (e.g. 



O i/H 2 0, N1/N2 or NH 3 ; see |Hollenbach fc Tielens] [1999) 



Inside of this photon-dominated layer, the grain temperatures are warm 
enough for CO to exist in the gaseous state (Tdust > 20 K), but water will 
remain as an ice. Thus, the C/O ratio ~ 1 leading to an active carbon-based 
chemistry in a vertical zone labeled as the "warm molecular layer" . The ion- 
molecule chemistry of this layer is powered predominantly by X-rays and 
Cosmic Rays (if present). Reactions between CO and molecular ions (pre- 
dominantly Hj) transfer a small fraction of this carbon into other simple and 



complex species (Aikawa et al. 19971. If CH4 is present on grain surfaces, 



and evaporates into the gas, it will also be a key pre-cursor to the creation of 
larger hydrocarbons and carbon chains. If Tdust < T su b for any product of this 
gas-phase chemistry, then that molecule freezes onto grains. In a large sense 
the gas-phase is acting as an engine for building complexity within molecular 
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Fig. 12. Chemical structure of protoplanetary disks (taken from Bergin et al. 2007 1. 
Vertically the disk is schematically divided into three zones: a photon-dominated 
layer, a warm molecular layer, and a midplane freeze-out layer. The CO freeze- 
out layer disappears at r < 30 — 60 AU as the mid-plane temperature increases 
inwards. Various non-thermal inputs, cosmic ray, UV, and X-ray drive chemical 
reactions. Viscous accretion and turbulence will transport the disk material both 
vertically and radially. The upper panels show the radial and vertical distribution 



of molecular abundances from a typical disk model at the midplane (Aikawa et al 
[1999} and r ~ 300 AU ||van Zadelhoff et al.||2003[). A sample of the hydrogen density 



and dust temperature at the same distance ( |D'Alessio et al!)|1999| ) is also provided 



ices (Tscharnuter & Gail 20071. As material advects inwards the region of 
the disk with Tdust > 20 K increases in depth and, in addition, the more 
tightly bound complex species created earlier (e.g. HCN, C2H2, C3H4) will 
eventually evaporate. An additional issue that is likely important is the cre- 
ation of frozen radicals via photodissociation on the grain surface. If the grain 
is warm enough these radicals can migrate and react on grain surfaces. This 



can lead to a large increase in molecular complexity (Garrod et al. 2008), 



perhaps contributing to organics detected in meteorites (see, e.g. Bernstein 
et al.[ 2002 and references therein) . 

Below the warm molecular layer lies the dense (n 3> 10 s cm -3 ), cold 



(T < 20 K) midplane where molecules are frozen on grain surfaces ( Willacy & 



L anger I 2000| Aikawa & Herbst 2001 1. This is the case for much of the outer 



disk (R ;> 20 - 40 AU) beyond the radial "snow-lines" . If ionizing agents are 
still present and T gas < 30 K, then the transition to the midplane and the 
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midplane itself are the main layers for deuterium fractionation in the disk 
( |Aikawa fe; Herbst" 20011. Gas phase deuterium fractionation requires heavy 
molecules to be present in the gas. In the limit of total heavy element (C, O, 
N) freeze-out, the chemi stry will be reduced to the fo llowing sequence Hg — > 



ff>D H 



Doff 



(Ceccarelli & Dominik 2005). 



Fig. 13 shows the CO abundance structure in a fiducial disk model to 



illustrate some of these effects more directly. The three panels refer to different 
models of dust settling, with lower values of the e parameter (defined in the 
figure caption) referring to a greater degree of dust settling. The top panel 
is for a disk with no settling. As can be seen the chemical structure divides 
into three layers as shown in Fig. 12 In addition the interior evaporation zone 
(inside the snow-line for CO) can be resolved. As grains settle into the disk 
the radiation penetrates more deeply as seen by the line delineating the r = 1 
surface and also in the chemical structure which moves the warm layer to 



lower depths (Jonkhcid et al. 2004) 



The case of water deserves special mention. Beyond the snow- line at ~ 1—5 



AU, water will exist mostly in the form of water ice (see Ciesla & Cuzzi 2006 



Davis 20071. Dust grains suspended in the atmosphere will typically not be 
heated above the evaporation temperature (Fig. [8]) and water will, for the 
most part, remain as ice. However, UV photons can photodesorb water ice 
from grain surfaces, producing a small layer where water exists with moderate 



Dominik et al. 2005). The 



abundance (xH 2 0,max ~ 10 , relative to ff; 
water vapor abundance peaks at the layer where the rate of photodesorption 
is balanced by photodissociation. The depth where the abundance reaches this 
maximum value depends on the strength of the UV field, the local density, 



and grain opacity (Hollcnbach et al. 2008 ). At lower depths photodissociation 



will destroy any desorbed ffO molecules and erode the mantle. At greater 
depths, the lack of UV photons leaves the ice mantle intact. Interior of snow 
line the water ice can evaporate from the grains and also be produced in the 
gas via the high temperature chemistry discussed in § |3.2| There is also a 
strong possibility of rapid movement of icy solids from the outer nebula that 



can seed additional water and other ices to inside the snow- line (Ciesla & 



Cuzzi 20061 



Radially there exist large gradients in the gas/solid ice ratio of molecular 
species in the midplane. The snow-line is species specific in the sense that CO 
should be present in the gas phase at greater radii than water (see Fig. [8]). 
This may not be a continuous transition with various species sublimating 
according to their molecular properties. Rather some molecules may be frozen 
in the water lattice (possibly clathrates in denser regions) and may evaporate 



along with water as seen in laboratory experiments ( Sandford & Allamandola 



1993 


Collings et al.| 2004 


Viti et al. 


2004) 



In the dense interior a wide range 
of high temperature kinetic reactions are likely active. These reactions can 



produce observed species, such as HCN and C2H2 (Gail 2002 Agundez et al. 



2008 1, and perhaps beyond. However, the eventual products are tempered 
by the destructive influence of photons. As noted in § |3.3[ thermodynamic 
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Fig. 13. Plot of the abundance of carbon monoxide as a function of radial and 
vertical height in three disk models with variable gas to dust ratios. The gas to dust 
ratio is parameterized by the parameter e, which is the gas to dust ratio in the upper 
layers relative to the gas to dust ratio in the midplane (see Calvet, this volume). 
Lower values of e are representative of the effects of the settling of dust grains. Also 
shown is the t — 1 surface for stellar photons which appears deeper in the disk as 
grains settle to the midplane. The top panel shows some of the basic disk chemical 
structure. Figure from Fogel et al. 2009, in preparation. 
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equilibrium can potentially be reached in the midplane of the inner few AU 
and in giant planet subnebula, with the removal of some elements from the gas 
as solids condense. However, the surface layers in the inner disk are dominated 
by the short timescale effects of stellar radiation. 

All of these processes must be viewed in light of the likelihood of dynam- 
ical mixing both radially and vertically. Models including mixing generally 
show that the vertical structure illustrated in Fig. 12 is preserved with some 

Beyond the in- 



widening of the warm molecular layer (Willacy ct al 



2006). 



ward advection of material the outward diffusion of hot gas from the inner 
nebula can bring gas to cold layers where molecules freeze onto grains surfaces 
( Tscharnuter & Gail 2007 ) . Moreover, the grains are evolving via coagulation 
and settling and the gas gradually becomes transparent to destructive radia- 
tion. Grain evolution ultimately results in the formation of planetesimals and 
planets (likely at different timescales for giant and terrestrial worlds) and the 
gas chemistry will continue until the gas disk dissipates. 



6 Summary 

We are approaching an age where studies of extra-solar systems can be in- 
formed by, and inform, studies of the chemical composition of bodies in our 
solar system. Some major gaps remain in our understanding that will benefit 
from this closer cooperation. These include (1) What is the relative impor- 
tance of thermodynamical equilibrium and kinetics in the inner disk; how 
does this inform the meteoritic inventory. (2) Inclusive of pre-solar grains, 
does any material, remain pristine and chemically unaltered from its origin 
in the parent molecular cloud? It is likely that the deuterium enrichments 
ultimately originate in the cold pre-stellar stages, but are these components 
dissociated and re-assembled into different form? (3) How deep inside the disk 
(radially) do cosmic rays penetrate and how might this influence dynamical 
evolution? (4) How extensive is the complex chemistry in the inner disk and 
how are these organics incorporated into planetesimals? (5) Diffusive mixing 
can be important, but how much and for how long is it active? These are just 
a sample of the many questions that remain. In this chapter we have outlined 
some of the basic physical and chemical processes that have laid a foundation 
for much of our current understanding. We hope that this chapter, and book, 
informs the next generation of researchers in order to untangle long-standing 
questions regarding the initial conditions, chemistry, and dynamics of planet 
formation, the origin of cometary ices, and, ultimately, a greater understand- 
ing of the organic content of gas/solid reservoirs that produced life at least 
once in the Galaxy. 
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